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1.0  INTRODUCTION  TO  THE  RESEARCH 

This  report  describes  a  continuation  of  research  into  the  modeling  of  water  particle 
freezing  for  application  to  adverse  weather  simulation  facilities.  The  research  was  initiated  in 
FY1996  to  investigate  the  physics  of  freezing  of  submillimeter  supercooled  water  particles  or 
droplets  in  both  natural  and  artificial  or  simulated  adverse  weather  environments.  The  first  phase 
of  the  research  was  reported  and  discussed  in  a  report  [1]  and  a  paper  [2].  The  work  has 
continued  into  FY1997  and  has  been  expanded  to  include  work  done  to  model  three-dimensional 
ice  accretions  on  surfaces,  as  well  as  modeling  the  near  field  of  water  spray  clouds  produced  by 
air-atomized  water  spray  nozzles.  Because  of  the  increased  scope  of  the  work,  a  single  report 
cannot  cover  all  of  the  work  phases.  Therefore,  the  present  report  covers  only  the  continued 
research  and  development  of  water  particulate  freezing  models  and  their  application  in  a  one¬ 
dimensional  multiphase  flow  code  to  predict  water  spray  freezing  in  ducted  air  flows. 

The  outline  of  the  report  is  as  follows: 

Section  2  provides  an  overview  of  applications  of  classical  homogeneous  nucleation 
theory  to  predict  supercooled  water  particle  freezing.  Specifically,  it  covers  and  discusses  the 
J(T)  and  I(T)  functions  (which  are  identical  except  for  nomenclature  used  by  various  authors) 
that  predict  the  rate-of-formation  of  freezing  centers  or  ice-germ  nuclei  that  initiated  the  freezing 
process. 

Section  3  describes  how  J(T)  and/or  I(T)  were  evaluated  in  the  present  study  and 
compared  to  their  evaluations  in  the  literature. 

Section  4  describes  how  J(T)  and/or  I(T)  were  modified  to  apply  to  the  heterogeneous 
freezing  of  nominal  purity  supercooled  water  droplets. 

Section  5  then  provides  calculated  results  from  incorporating  the  modified  heterogeneous 
droplet  freezing  model  into  the  AEDC  1DMP  code  and  evaluating  some  of  the  same  flow  cases 
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as  described  in  references  [1,2]  for  a  ten-droplet  water  spray.  The  results  from  a  previously 
developed  simple  heterogeneous  freezing  model  and  the  new  modified  homogeneous  nucleation 
freezing  function  model  are  compared  and  evaluated. 

Section  6  then  summarizes  the  current  status  of  research  and  development  under  this 
phase  of  the  research  task  investigation,  provides  conclusions  about  accomplished  and  needed 
future  research,  and  makes  recommendations  about  how  to  pursue  future  development  of  adverse 
weather  simulation  test  capabilities  in  ground  test  facilities. 

2.0  OVERVIEW  OF  CLASSICAL  HOMOGENEOUS  NUCLEATION  THEORY 

2.1  Nucleation  Probability  Functions 

This  section  very  briefly  reviews  the  homogeneous  nucleation  probability  functions  of 
the  literature  and  describes  the  selection  process  used  to  select  one  of  the  functions  for 
application  to  development  of  a  heterogeneous  freezing  model  for  supercooled  water  droplets. 

In  this  context,  homogeneous  nucleation  refers  to  the  creation  of  a  condensed  phase  of  a 
pure  substance  by  the  thermodynamically  controlled  formation  of  condensed  phase  precursors. 
These  precursors  are  metastable  clusters  of  molecules  that  potentially  can  serve  to  trigger  the 
formation  of  the  condensed  phase,  that  is,  act  as  condensation  sites.  Specifically,  the  precursors 
can  be  thought  of  as  predroplets,  or  nandroplets,  when  a  liquid  condenses  from  a  supersaturated 
vapor,  or  as  ice  germs  when  ice  crystals  form  within  a  liquid.  Heterogeneous  nucleation  follows 
the  same  thermodynamic  process  of  molecular  cluster  formation  as  occurs  in  homogeneous 
nucleation,  except  that  some  foreign  material,  such  as  a  particle  or  a  surface,  reduces  the  energy 
required  to  form  the  precursors.  Thus,  heterogeneous  nucleation  occurs  at  higher  phase 
temperatures,  in  general,  than  homogeneous  nucleation. 
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The  literature  on  nucleation  and  condensation  processes,  that  start  with  the  formation  of 
metastable  molecular  clusters  in  the  non-condensed  phase  of  a  substance,  is  reasonably  large  in 
both  papers  [3-38],  and  books  [39-43].  The  “condensed”  phase  may  be  a  liquid  droplet 
condensing  from  a  (non-condensed)  supersaturated  vapor,  or  it  may  represent  a  proto-crystal  or 
ice  germ,  in  the  case  where  pure  supercooled  liquid  water  is  beginning  its  spontaneous  freezing 
process,  at  some  temperature  below  273.16  K.  The  freezing  process  that  begins  to  occur  in 
sprays  of  supercooled  liquid  metal  droplets  is  essentially  the  same  process  as  that  in  which  pure 
supercooled  water  spray  droplets  freeze,  except  that  liquid  metals  can  be  alloys  and  can  freeze  in 
different  solid  phases,  depending  on  drop  size,  composition,  and  temperature  at  the  time  of 
solidification  nucleation.  However,  pure  liquid  water  droplets,  supercooled  and  freezing,  turn 
into  ice  I-type  solid  crystalline  particles  under  normal  “atmospheric”  conditions. 

As  stated,  the  purpose  of  this  section  of  the  report  is  to  briefly  overview  the  probability 
functions,  denoted  in  the  literature  as  J(T)  or  I(T),  which  predict  the  production  rate  of  nucleation 
sites  per  unit  volume  per  unit  time.  In  the  present  research  effort,  the  theory  for  homogeneous 
nucleation  of  “pure”  liquid  water  droplets,  or  homogeneous  freezing,  has  been  modified  to 
account  for  heterogeneous  (or  mote-induced)  freezing  of  water  droplets.  Heterogeneous  freezing 
is  the  usual,  typical,  or  most  probable  freezing  process  for  liquid  water  of  average  purity,  such  as 
either  unfiltered,  or  even  filtered,  and  distilled  drinking  water.  The  modified  homogeneous 
nucleation  freezing  theory  can  be  applied  in  a  straight-forward  way  to  general  numerical  spray 
computation  codes  (numerical  spray  CFD  codes)  that  track  the  position,  lifetime,  kinetic,  and 
thermal  states  of  particles  or  particle  group  packets  that  are  injected  and  convected  in  gas  flows. 
Moreover,  the  modified  theory  can  incorporate  probabilistic  functions  that  extend  the  theory  to 
stochastic  processes  of  freezing.  To  explain  the  modified  theory,  it  is  necessary  to  review  or 
overview  classical  homogeneous  nucleation  theory  for  particulate  freezing.  The  theory  begins  in 


Germany  with  the  work  of  Volmer,  et  al.  [3],  Becker,  et  al.  [4],  and  others  with  further 
development  by  Turnbull,  and  co-workers,  e.g.  [7,8].  Its  exposition,  however,  is  begun  with 
Mason  [14],  for  convenience. 


2.2  Mason  (1952) 
In  his  paper,  Mason 


[14]  provides  J(T)  in  the  form 


(1) 


Mason  then  evaluates  log  J(T) 1  as 


log  7(r)  =  32.84  +  log  T- 


U 


760c4 


2.30  kT  {Tb-T)2T 


(2) 


In  these  expressions,  the  parameters  are: 

n  =  number  of  molecules  per  cubic  centimeter  of  the  condensed  phase, 
k  =  Boltzmann’s  constant, 
h  =  Planck’s  constant, 

T  =  the  absolute  temperature,  Kelvins, 

Ta  =  the  bulk  water  freezing  temperature,  Kelvins  (273.1 6K) 

U  s  the  “activation  energy  for  self-diffusion  of  a  molecule  in  the”  non-condensed  phase, 
Wc  =  “the  work  of  nucleus  formation”. 

Mason  provides  Wc  as 

Wc=^(TSLA  =  ^cSLQ)rc2  (3) 

attributing  this  to  Frenkel  [39].  These  parameters  are: 


1  Footnote:  In  ( )  refers  to  the  natural  log,  whereas  log  ( )  refers  to  log  base  10. 
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oSL  =  the  “specific  surface  energy  of  the  crystal-liquid  interface”. 

A  —  the  surface  area  of  the  nucleus  (the  molecular  cluster,  or  ice  germ  precursor). 
co  m  shape  factor  for  nucleus,  relating  the  surface  to  volume  ratio  for  the  shape  of  the 
cluster  (co  is  approximately  21  -  23  in  value), 
r.  =  the  radius  of  a  sphere  inscribed  in  a  nucleus  of  critical  size,  i.e.,  that  sized  nucleus 
that  can  form  a  nucleation  or  crystallization  site. 

Mason  evaluates  co  by  assuming  that  the  nucleus  forms  a  hexagonal  prism  with  height  equal  to 


twice  the  radius  to  a  prismatic  side  (see  sketch  below): 

Schematic  of  Prismatic  Ice  Germ  Nucleus 


Top  or  Bottom  View  Side  View 


The  volume  of  this  hexagonal  prism  is  V  =  4^3  rc3  =  6.9282rc3 

(4) 

The  surface  area  is  A  =  12^/3~ rc2  =  G)rc2  =  20.785r2 

(5) 

.  .  A  12V3  r]  3 

The  surface  to  volume  ratio  is  =  r-  ,  - 

V  4V3  rc3  rc 

(6) 

V  V  V 

Note  also  that  A  can  be  written  in  terms  of  V  as  A  =  3/6  4/3  V/3 

(7)' 

So  that  the  derivative  dA/dV  is 
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(8) 


f-A  -.3%  /|K  y"/4  _  2-3^  .4^ 

dV  1 3  J 

Substituting  for  V 

^■  =  2-3^-4^/(4V3rIsf 


(9) 


dA  _  2  3^  _  2 
"  3%  rc  "  rc 


(10) 


This  parameter,  dA/dV,  is  important  because  Mason,  in  a  later  paper  [25(1960)],  shows  that  the 
change  in  Gibbs  Free  Energy  for  a  collection  containing  g  molecules  that  have  passed  from  the 
liquid  phase  to  form  a  cluster  in  the  solid  phase  is 

AG  =  (/is  - [iL)g  +  A<Jsl  (11) 


where 

fiL  =  Gibbs  Free  Energy  of  molecule  in  liquid  phase, 

Hs  =  Gibbs  Free  Energy  of  molecule  in  an  ice-like  cluster, 

g  =  number  of  molecules  in  the  cluster. 

Mason  [25(1960)]  defines  the  critical  cluster  of  molecules  as  the  metastable  cluster  reached 
when  AG  achieves  a  maximum,  thus  the  two  conditions 


dg 


(12) 


and 

AG  =  AGmax  (13) 

define  the  cluster  size  when  a  metastable  ice  germ  (cluster)  has  formed. 

Therefore,  for  this  condition, 
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(14) 


£ 


But, 


and 


dAG 

dg 


=  0  =  {^s~Hl)+ 


dA 


dA  =  dA  dV_ 
dg  dV  dg 


(15) 


V  =  Xg 


(16) 


where  X  =  the  effective  volume,  per  molecule,  in  the  cluster. 

Thus, 

(H) 

dg  dV 


An  evaluation  of  X  gives 


A  _ 

NAPs 


(18) 


where 

=  the  molecular  weight  of  ice, 

Na  s  Avogadro’s  number, 

ps  =  the  density  of  ice. 

Thus, 

dA  _dA  Wmol 
dg  dV  NAp,  ' 


Consequently,  the  condition  for  metastability  of  die  cluster  is 


W~Ps)  = 


OsF-i  dA 
NaPs  dV' 


(19) 


(20) 
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However,  Mason,  e.g.  [1960],  shows  that 


where 

L*f  =  the  heat  of  fusion,  per  molecule  (treated  as  a  negative  quantity), 

SL  =  the  entropy  of  a  molecule  in  the  liquid  phase, 

S  =  the  entropy  of  a  molecule  in  the  ice-like  cluster  or  solid  phase. 

Thus,  by  integration 

T„  1 

and,  by  substitution,  therefore,  based  on  equation  (20), 

y _<^mot  dA_\llTdT 

L  *s)  NAp,  dv  It’ 


Substituting  for  dA/dV 

Qsfflm*  2  _\L\dT 
NaPs  rc  l  T 

Solving  for  rc,  and  rearranging  terms  gives  rc  as 


But,  the  heat  of  fusion  per  unit  mass  of  the  liquid  is 


Therefore,  the  equation  for  the  critical  radius  of  the  metastable  molecular  cluster  is 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 
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(27) 


Therefore,  Wc  can  be  written  as 


16 So], 


(31) 


However,  in  his  1952  paper  [14],  Mason  simply  introduces  an  equation  for  rc  that  Mason  calls 
“Thompson’s  equation,”  without  reference,  as 


^  &SL 

PA 


Jo 

T0-T' 


(32) 


Mason  later  provides  a  derivation  of  this  equation  in  his  1957  text,  The  Physics  of  Clouds,  [44]. 
Clearly,  equation  (32)  is  also  obtainable  from  equation  (27)  when  LF  is  treated  as  a  constant,  by 
expanding  and  linearizing  the  resulting  log  function.  In  the  previous  equation,  ps  is  the  density 
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of  the  condensed  phase,  here  assumed  to  be  Ice-I,  normal  ice.  By  using  equations  (3),  (30),  and 
(32)  in  equation  (1),  we  can  write  the  equation  for  J(T)  as 


nkT  U  \y3P*.a(  2gsJ(  T0 

r.-r 


By  using  T0  =  273.16  K,  and  simplifying,  we  get 


nkT  [  U  4  gsl3co  (  273.16 


T  (rf,  \  W  -r  ''SL  |  -  - 

j(r)____exp  ^  3kTP2sL2F{213A6-T 


Equation  (2)  is  equation  (34),  after  some  of  the  parameters  have  been  numerically 
evaluated.  Mason  provides  values,  in  his  1952  paper,  for  the  following  parameters  in  J(T),  as 


follows. 


The  interfacial  surface  tension,  GSL,  was  evaluated  by  Mason  as  the  difference  between 
the  surface  tension  of  ice  at  -40°C  and  of  liquid  water  at  -40  C,  that  is, 

102^-80^ 

SL  cm  cm 


^=22-^  =  0.022^-. 

SL  cm2  m 

Gs  was  estimated  from  breaking  hydrogen  bonds  normal  to  a  selected  face  of  an  ice  crystal. 

The  value  for  U  was  estimated  at  U  =  3.3 x  10'13 erg,  or  3.3x10  7,for0<T<-10  C, 

and  this  same  value  taken  at  T  =  -40°C.  The  density  of  ice  and  its  latent  heat  of  fusion  were 
taken  as 

ps  -  0.92  gm/cm2  =  920  kg/m3 
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Lf  =  -3.33  xl09H£  =  -3.33  *  105  — 
gm  kg 

The  shape  factor,  co,  was  taken  by  Mason  to  be  (0  =  23.  (However,  it  was  computed  in 
the  present  study  as  about  21  instead.) 

The  next  parameter  is  the  number  n  of  molecules  of  liquid  water  per  unit  volume.  This 
p  -N 

parameter  is  given  by  n  =  — — —  where  pL  is  the  liquid  density,  NA  is  Avogadro’s  number,  and 


W 


mol 


is  the  molecular  weight  of  water.  Thus,  for 


p,= i-n 

cm 


Na  =  6.022169  x  10 


23  molecules 


gm  -  mole 
Wmol  =  1 8  gm/gm  -  mole, 

n  is  computed  as  n  =  3.3456  x  10  - ^ —  =  3.3456  x  10  - r - . 

cm  m 

The  remaining  parameters  are  Boltzmann’s  constant,  k,  and  Planck’s  constant,  h,  given  by 

k  =  1.380622  x  10’16  —  =  1.380622  x  10"23  — 

K  K 

and 

h  =  6.626196  x  10_27erg  •  s  =  6.626196  x  10-347  •  j  . 


Evaluation  of  log  J(T) 

Using  these  values  for  the  parameters  in  the  J(T)  equation  Mason  evaluated  log  J(T)  as 
log  J  (T)  =  log(3.3456  x  1022)  +  log(l. 380622  x  10'16)  +  log!  -log(6.626196  x  10"27) 


U 


4  1 


cr3  (23)(273.16)2 


1 


(35) 


2.303  kT  3  2.303  (1 . 380622  xlO-16)  T  (0.92)2 (3.33  xlO9)2  ( T0-T ) 
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Evaluating  this,  we  get 


log  7  (1)  =  32.843  +  log  T- 


U  766.8(7 ]L 
2.303kT  T(T0  -  T)2 


(36) 


where  cSL  is  in  erg/cm2  units  and  J(T)  would  be  in  units  of  nuclei  per  cubic  centimeter  per 
second.  With  minor  numerical  differences,  this  is  Mason’s  equation,  equation  (2). 

In  the  present  study,  using  the  values  of  the  parameters  as  given  in  SI  units,  and  using 
natural  logarithms,  this  equation  has  been  put  into  the  form 


In  7(7)  =  89.440  + In  7- 
Altemately,  this  equation  gives  J(T)  as 

J(T)  =  6.9708  x  1038r  exp 


2390.2 

252  f 

7 

T  { 

2390.2 

252  f 

T 

T  { 

273.16-7 


273.16-7 


J  J 


(37) 


(38) 


with  units  of  nuclei  per  meter  cubed  per  second. 

In  the  following  subsections  of  the  report,  alternate  forms  for  J(T)  are  described  which 
are  given  in  the  literature.  J(T)  has  been  also  given  the  symbol  I(T).  In  this  report,  J(T)  and  I(T) 
are  exactly  the  same  functions. 

2.3  McDonald  (19531 

In  this  paper,  McDonald  [18]  provides  a  critical  review  of  homogeneous  nucleation 
theory.  In  particular,  with  reference  to  the  application  of  the  theory  to  predict  the  spontaneous 
freezing  of  supercooled  submillimeter  water  droplets  by  homogeneous  nucleation,  McDonald 
reviews  and  assesses  the  values  of  the  physical  parameters  in  the  argument  of  the  function,  J(T). 
McDonald  writes  the  7(7)  function  as  (with  a  sign  correction). 
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nkT 


{_ 

(  A  +  Fc') 

1 

kT 

\  K1  ) 

(39) 


A  is  the  activation  energy  for  self-diffusion  of  liquid  molecules  near  the  water-ice  germ  interface. 
Fc  is  the  work  of  formation  of  the  ice  germ  or  molecular  cluster,  n  is  the  number  of  molecules 
per  unit  volume  in  the  liquid  phase,  k  is  Boltzmann’s  constant,  and  h  is  Planck’s  constant. 
McDonald  gives  Fc  as 

Fc=osgr?/3  (40) 

where  rc  is  the  critical  radius  of  the  molecular  cluster,  given  by 


rc  = 


2oT 

s  o 


(41) 


In  these  equations,  ex  is  the  specific  surface  free  energy  of  the  water-ice  germ  interface,  g  is  a 
geometric  factor  such  that  g  rc2  is  the  total  surface  area  of  the  ice  germ  or  critical  embryo,  ps  is 


the  density  of  ice,  LF  is  the  latent  heat  of  fusion  of  ice,  and  To  is  the  melting  temperature  of  ice, 
273.16  K.  Thus,  the  expression  for  J(T)  is  essentially  the  same  as  that  presented  earlier  by 
Mason  and  co-workers. 

The  main  thrust  of  McDonald’s  investigation  is  an  examination  of  the'  parameters  A,  <JS, 
and  Lr  McDonald  reviews  the  work  of  others,  including  Mason  [14],  in  an  attempt  to  obtain 
thermodynamically  correct  values  of  these  parameters,  in  terms  of  the  best  theoretical  model  of 
an  ice  germ  structure.  McDonald  notes  that  Lr  the  heat  of  fusion  of  ice  is  not  constant,  but 
decreases  with  decreasing  temperature.  (This  was  discussed  by  the  present  author  in  previous 
work  [1,2].)  McDonald  also  finds  both  A  and  as  to  vary  with  temperature. 

McDonald’s  rigorous  analysis  is  well  worth  examination  by  those  with  expertise  in 
molecular  physics  because  the  levels  of  uncertainty  in  the  estimations  or  calculations  of  A,  <JS, 
and  Lr  that  McDonald  demonstrates,  implies  that  serious  further  theoretical  and  experimental 
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research  is  required  to  put  homogeneous  nucleation  theory  on  a  rigorous  basis.  This  is  one  of 
McDonald’s  main  conclusions  in  1953  and  it  appears  to  be  valid  even  at  present.  In  actual 
practice,  therefore,  the  application  of  the  freezing  function  J(T),  by  Mason  and  others,  is  a 
correlative  theory,  rather  than  a  predictive  theory,  due  to  the  uncertainties  in  the  physical 
parameters  of  the  theory. 

2.4  Dav  0958} 


J.  A.  Day  [21]  presented  a  functional  form  for  J(T)  that  he  attributed  to  J.  B.  McDonald 
[18].  This  equation  is 


r,_.  nkT 
7(7)  =  — -exp 
h 


30.7(7 


SL 


kT 


kp)L\T 


f 


In  — 
T 


(42) 


where  T0  and  the  other  parameters  have  been  defined  in  the  previous  section.  Obviously,  this 
equation  can  be  developed  from  equation  (1)  with  Wc  evaluated  from  equation  (31)  and  with  Lf 
treated  as  constant. 

2.5  Mason  (1958) 

Mason  presented  another  version  of  the  J(T)  equation  in  this  1958  paper  [22].  In  this 
paper,  Mason  uses  I(T)  to  represent  the  nuclei  density  rate  equation.  The  equation  provided  was 


U  l.llxlO17^ 


log  I(T)  =  32.84  +  log  T  - 


2.303 kT 


L\T 


( 


V  ’ 


(43) 


In  — 
T 


\  J 

The  same  comments  apply  here  to  this  equation  as  for  the  equation  Day  presented. 
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2.6  Lansham  and  Mason  (1958) 

In  1958,  Langham  and  Mason  reviewed  the  heterogeneous  and  homogeneous  nucleation 
of  supercooled  water  [23]  and  presented  a  functional  equation  for  1(T)  that  is  the  same  as  that  of 
Mason  (1958)  [22]. 

2.1  Mason  (I960) 

In  a  1960  paper,  Mason  [25]  applied  nucleation  theory  to  predict  both  the  formation  of 
water  aerosol  droplets  from  supersaturated  water  vapor,  and  also,  homogeneous  nucleation  or 
freezing  of  supercooled  water  droplets.  Mason  also  reviews  the  formation  of  water  droplets  by 
condensation  of  vapor  on  foreign  nuclei  such  as  ions,  hygroscopic  and  non-hygroscopic  motes, 
and  on  “mixed”  nuclei. 

Mason  attributes  the  theoretical  development  of  condensation  theory  to  Becker  and 
Doring  [4],  Zeldovich  [45],  and  Turnbull  and  Fisher  [7].  Mason  discusses  the  validity  of  the 
theory,  as  well  as  the  validity  of  the  experiments  that  provided  data  on  droplet  condensation  and 
freezing. 

Mason  discusses  the  freezing  of  submillimeter  supercooled  droplets  by  both 
heterogeneous  and  homogeneous  nucleation  processes  and  the  applicable  theories  for  both. 
Mason  also  reviews  the  effectiveness  or  activity  of  various  mineral  substances  to  act  as  ice 
nucleating  motes.  The  most  active  substance  Mason  found  was  silver  iodide  which  Mason  says 
produces  one  ice  crystal  per  10000  Agl  particles.  This  statement  has  relevant  implications, 
therefore,  in  the  interpretation  of  J(T).  It  implies  that  the  population  number  densities  for 
nucleation,  as  correlated  by  J(T)  (or  I(T))  must  reflect  only  the  active  or  effective  nuclei,  or  the 
active/effective  motes,  and  not  the  actual  or  physical  nuclei/mote  population  number  densities  or 
their  rates  of  creation.  Especially  this  will  be  true  if  J(T)  is  used  to  predict  heterogeneous 
freezing  by  correlating  the  J(T)  function  (actually  csl)  to  a  given  set  of  experimental  data. 
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A  derivation  of  1(T)  was  provided  by  Mason  in  this  paper  [25]  that  leads  to  the  equation 


log  7(1)  =  32.84  +  log  T- 


U 


16^3  a 


SL 


2.303 kT 


2.303 p]  kT 


£LtdT 


-t2 


(44) 


It  is  clear  that  this  equation  can  be  derived  from  the  approach  described  in  section  2.2  of 
the  present  study.  In  discussing  this  equation,  Mason  notes  that  <JSL  probably  cannot  be  calculated 
with  sufficient  accuracy  to  permit  using  I(T)  to  predict  freezing  nucleation  rates.  However,  he 
also  notes,  that  by  using  experimental  data  for  the  freezing  of  pure  supercooled  water  droplets, 
oSL  values  can  be  obtained  by  correlating  I(T)  with  appropriate  sized  drops  freezing  at  measured 
temperatures.  If  the  presence  of  motes  in  ordinary  water  acts  to  change  the  crSL  values,  then  it 
should  also  be  possible  to  correlate  I(T)  to  the  heterogeneous  freezing  of  supercooled  droplets 
and  obtain  gsl  versus  drop  size  for  each  set  of  experimental  data.  This  has  been  done  in  the 
present  research  investigation  to  obtain  a  J(T)  function  applicable  to  heterogeneous  nucleation 
freezing  of  supercooled  water  droplets  and  will  be  presented  in  section  4. 

2.8  Fletcher  (196(7) 

In  this  paper  [24],  Fletcher  writes  the  nucleation  site  population  density  function  as 

J(T)  =  £exp{-  AG*/kT}  (45) 

where  “ K  is  a  kinetic  constant  typically  of  order  10  cm  sec  for  the  cases  which  we  shall 
consider,”  k  is  Boltzmann’s  constant,  and  T  is  absolute  temperature,  as  utilized  previously. 
Fletcher  discusses  the  extension  or  application  of  nucleation/condensation  theory  to  account  for 
the  effect  of  foreign  surfaces  or  particles  in  modifying  (reducing)  the  value  of  the  work  of 
formation  of  the  molecular  cluster  or  embryo  nuclei,  denoted  AG*.  Thus,  if  a  theoretical  model 
can  be  developed  of  how  substances,  either  as  motes  or  surfaces,  can  change  the  value  of  AG*, 
then  this  model  can  be  incorporated  into  J(T)  to  create  a  heterogeneous  nucleation  site  population 
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rate  equation.  Theoretically,  therefore,  as  Becker,  Turnbull,  Mason,  and  others  have  also 
discussed,  heterogeneous  nucleation  freezing  of  supercooled  water  droplets  can  be  predicted  with 
the  same  theoretical  tools  used  to  predict  homogeneous  freezing  of  submillimeter  supercooled 
water  droplets. 

Fletcher  discusses  the  effects  of  foreign  surface  geometry  and  size  on  AG*  for 
heterogeneous  nucleation  or  the  formation  of  the  ice  germ  or  embryo.  He  presents  a  function, 
f(m,x),  that  modifies  the  free  energy  of  formation  for  homogeneous  nucleation,  i.e.,  Fletcher 
gives 

AG*  =  AG'0  f(m,x)  (46) 

where  AG*  is  the  free  energy  of  formation  of  a  nucleus  of  critical  radius  under  homogeneous 

nucleation  conditions,  such  as  given  by  equation  (31)  and  where 

m  =  cos  9 .  (47) 

Cosine  6  is  the  cosine  of  the  contact  angle  between  ice  germ  (embryo)  and  the  foreign  particle 
and  x  is  defined  as 

x  =  R/r*.  (48) 

R  is  the  radius  of  a  spherical  foreign  particle  and  r*  is  the  radius  of  curvature  of  the  surface  of  the 
critical  embryo. 

Fletcher  [24]  provides  distributions  of  f(m,x)  in  a  figure  (Fig.  1  of  [24])  obtained  from 
solving  the  “nucleation  equations  for  a  spherical  cap  embryo  nucleating  upon  a  perfect  spherical 
particle.”  Details  of  this  model  are  provided  by  Fletcher  in  another  paper. 

Fletcher  shows  that  this  theoretical  model  for  accounting  for  the  effects  of  foreign 
particles  on  the  freezing  of  supercooled  water  predicts  the  correct  trends  when  applied  to  silver 
iodide  particles  inducing  higher  spontaneous  freezing  temperatures  in  supercooled  droplets  of 
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water.  Thus,  the  use  of  modified  homogeneous  freezing  theory  to  predict  heterogeneous  freezing 
can  be  justified  by  more  than  one  theoretical  approach. 

2.9  Kuhns  and  Mason  Cl 968) 

This  paper  [29]  continues  the  work  of  Mason  on  the  assessment  and  theoretical 
evaluation  of  the  freezing  of  submillimeter  supercooled  droplets  of  pure  water.  Kuhns  and 
Mason  review  and  describe  the  work  of  Mason  (1957,  1958,  1960)  and  others  who  have  studied 
the  freezing  of  small  water  droplets.  Kuhns  and  Mason  report  experimental  data  they  obtained 
on  freezing  of  freely  falling  water  drops,  noting  that  much  of  the  existing  data  to  1968  consisted 
of  freezing  droplets  that  were  supported  or  in  contact  with  surfaces,  wires,  hypodermic  needles, 
etc.  Much  of  this  paper  describes  the  apparatus  used  by  Kuhns  and  Mason  to  perform  their 
experiments  and  to  analyze  the  freely  falling  droplets  to  estimate  their  temperature  time  history, 
as  well  as,  size.  They  also  investigated  the  effect  of  different  ambient  gases  surrounding  the 
particles  on  the  particle  freezing  process. 

These  authors  review,  and  apparently  expand,  the  heterogeneous  freezing  theory  of  Bigg 
[16],  as  well  as,  the  homogeneous  theory  developed  by  Mason  (1952,  1958,  1960)  from  the 
Turnbull-Fisher  work  [7].  They  report  the  same  J(T)  function  that  Mason  presented  in  his  1960 
paper  (see  section  2.7  of  the  present  report).  Based  on  the  probability  theoiy  of  Bigg  [16],  and 
using  the  concept  of  a  “median”  freezing  temperature  for  a  group  of  drops  all  of  the  same  size, 
but  freezing  stochastically  about  this  median  temperature,  the  authors  define  the  “median” 
freezing  event  as  occurring  when  the  product  of  J(TF),  particle  volume,  V,  and  time,  t,  reach  a 

critical  value.  That  is,  when  the  freezing  probability,  Pr  reaches  PF  =  — ,  half  of  the  drops  are 

2 

assumed  frozen,  [16].  To  illustrate  this  theory,  let  the  probability  of  a  drop  remaining  liquid  be 
given  by 
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The  present  author  utilizes  equation  (49)  differently.  Equation  (49)  can  be  interpreted  or 
utilized  in  several  ways.  With  aSL  determined,  equation  (49)  can  be  used  to  predict  the  freezing 
temperature,  TF,  of  different  sized  droplets.  Alternately,  if  the  “median”  freezing  temperature, 
Tf,  of  a  given  sample  of  drops  of  all  the  same  size  is  known,  then  equation  (49)  can  be  used  to 
evaluate  a  “median”  value  for  aSL.  However,  the  present  author  interprets  equation  (49) 
somewhat  differently. 

There  are  two  reasons  for  the  present  author  using  a  different  interpretation  to  the 
product,  JVt.  From  a  physical  interpretation,  the  number,  Nc,  of  (active)  freezing  nuclei  present 
in  a  single  sample  or  droplet  of  supercooled  water  at  time,  t,  must  be  given  by  the  following 
integral 

Nc{t)=]j{T,D,t)v{D,t)dt  (50) 

* o 

where  tF  is  the  time  at  which  the  particle  begins  to  freeze,  and  t0  is  the  time  at  which  the  particle 
first  began  to  supercool  (the  time  when  its  temperature  reached  273. 16K).  In  this  integral,  the 
function,  J,  which  gives  the  population  net  rate-of-creation,  has  been  generalized  to  be  a  function 
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not  only  of  particle  temperature,  T,  but  also  of  particle  size,  given  by  its  diameter,  D,  and  time,  t. 
Thus,  equation  (50)  can  apply  equally  well  to  heterogeneous  nucleation  events  in  supercooled 
water. 

With  this  interpretation,  the  minimum  or  critical  value  for  Nc  is  unity,  that  is,  for  a  single 
drop  undergoing  freezing, 

(51) 

Hence,  by  the  present  author’s  interpretation  at  least  one  active  freezing  site  must  come  into 
existence  in  each  sample  to  initiate  freezing. 

The  second  reason  for  this  interpretation  is  that  the  freezing  model  must  be  implemented, 
in  the  present  study,  into  a  code  wherein  all  droplets  of  the  same  size  freeze  at  the  same  time.  At 
present,  the  limitations  of  the  code  preclude  modeling  wherein  half  of  the  droplets  of  a  given  size 
freeze  at  the  median  freezing  temperature  and  half  remain  liquid.  Hence,  when  a  defined 
measure  of  the  product  JVt  reaches  a  critical  value  for  the  droplets  of  that  size  class,  assumed  in 
this  study  to  be  unity,  all  of  these  droplets  begin  to  freeze.  The  measure  of  JVt  used  in  the 
present  study  is  an  integral  of  the  product  JV  dt,  as  defined  below. 

In  the  present  study,  a  criterion  or  a  measure  of  JVt  for  both  homogeneous  and 
heterogeneous  nucleation  freezing  was  adopted,  given  by  the  integral  equality 

]j{T,D,t)V(D,t)dt  =  l  .0  (52) 

When  this  integral  reaches  unity,  all  drops  of  size  D  are  assumed  to  begin  freezing.  This 
criterion  (equation  (52))  was  adopted  and  has  been  applied  in  the  present  research  study  to 
predict  heterogeneous  nucleation  freezing  of  water  spray  droplets  in  ducted  air  flow,  by  being 
incorporated  into  a  numerical  model  for  one-dimensional,  multiphase  flow  described  in  a 
previous  report  [1],  and  paper  [2].  Details  of  this  model  are  presented  later  in  the  present  report. 
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Returning  to  a  discussion  of  the  Kuhns  and  Mason  paper,  Kuhns  and  Mason  determine 
the  value  of  aSL  required  in  J(T)  by  requiring  20-micron  drops  to  freeze  at  a  median  freezing 


temperature  of  -37C  after  a  period  of  1  second.  This  resulted  in  a  cSL  value  of 

aSL  =  19.7 =  0.0197  V  2  (53) 

cm  / m 

where  oSL  is  the  same  free  energy  of  the  ice-water  interface  or  surface  as  described  previously. 
With  this  set  of  conditions,  with  aSL  of  19.7  erg/cm2,  the  other  freezing  temperatures  of  droplets 
of  different  sizes  can  be  determined  when 

7(7>()V;.(lsec)  =  0.7  (54) 


where 


spontaneous  (median)  freezing  temperature  of  droplets  of  diameter,  di 
7C 

—df,  the  droplet  volume  for  a  drop  of  diameter,  dr 


They  presented  the  resulting  TF_  vs  d.  curve  and  compared  experimental  data  from  various  sources 

against  it,  including  their  own  data.  The  results  of  the  comparison  were  good. 

Kuhns  and  Mason  also  estimate  the  size  of  the  critical  embryo,  or  nucleus  or  ice  germ 
that  reaches  a  metestable  equilibrium,  and  thus  has  an  opportunity  to  nucleate  the  freezing 
process.  The  critical  radius  of  the  cluster  was  given  by  the  same  equation  as  equation  (27), 
namely 


where,  as  before, 

ps  =  the  density  of  ice, 


2(7 


rc  = 


SL 


P'\ 


r-  Lf  dT 


(55) 
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L  m  the  heat  of  fusion  (or  melting)  of  ice,  (now,  however,  treated  as  positive). 


If  Lf  is  constant,  rc  is  given  by 


2cr 


rc=- 


SL 


PsLMf- 


(56) 


By  assuming  different  values  of  oSL,  the  corresponding  values  of  the  critical  radius  rc  were 
computed,  and  the  number  of  water  molecules  that  form  the  ice  precursor  or  embryos  were 
estimated.  The  authors  obtained  150-300  water  molecules  depending  on  cSL  in  the  temperature 
range  investigated.  The  authors  also  showed,  albeit  with  some  approximation,  that  the  results  of 
a  statistical  thermodynamics  model  of  water  molecule  aggregates,  termed  the  “flickering  cluster 
model,  explored  by  Nemethy  and  Scheraga  [46],  can  be  interpreted  to  predict  the  same  number 
of  molecules  required  for  a  meta-molecular  cluster  and  for  a  Gsl  value  of  about  20  erg/cm  .  Thus, 
it  is  implied  that  the  methods  of  statistical  thermodynamics,  applied  to  a  special  model  of 
molecular  aggregates,  could  be  used  to  predict  the  spontaneous  freezing  of  pure  supercooled 


water. 


2.10  Anderson.  Miller.  Kassner,  Jr.,  and  Hagen  (1980) 

This  paper  [33]  reviews  condensation-freezing  nucleation  of  small  water  droplets  in  an 
expansion  cloud  chamber.  The  authors  report  on  their  experimental  cloud  chamber  test  results 
where  submillimeter  liquid  droplets  first  nucleate  from  a  supersaturated  vapor  phase  then 
undergo  spontaneous  homogeneous  (nucleation)  freezing.  The  freezing  process  is  reported  to 
occur  where  chamber  temperature  is  “near  -40°C.”  The  authors  also  note  minor  effects  of  an 
electric  field  applied  to  the  cloud  chamber  which  reduces  the  presence  of  ions  on  which  droplets 
can  nucleate  and  thus  freeze. 
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The  authors  investigate  whether  the  occurrence  of  frozen  nuclei  in  cold  supersaturated 
vapors  is  due  to  the  direct  formation  of  ice  germs  or  ice  nuclei  by  vapor  to  solid  nucleation 
processes  or  whether  it  occurs  due  to  a  two-step  process  of  nucleation  of  vapor  first  to  form 
liquid  droplets,  then,  second,  for  the  liquid  droplets  to  freeze  by  liquid  to  solid  nucleation.  An 
understanding  of  which  of  these  processes  can  account  for  the  presence  of  atmospheric  ice 
particles  would  improve  our  general  knowledge  about  the  mechanisms  of  formation  of  arctic 
precipitation  clouds. 

The  authors  discuss  the  literature  on  water  nucleation,  including  vapor  to  liquid,  vapor  to 
solid,  and  liquid  to  solid  nucleation.  They  provide  expressions  for  two  nucleation  site  creation 
rate  functions,  JJT)  and  JJT).  JVL(T)  is  the  homogeneous  nucleation  site  creation  rate  of  liquid 
water  droplets  forming  by  nucleation  from  supersaturated  water  vapor.  JVS(T)  is  the 
homogeneous  nucleation  site  creation  rate,  per  unit  volume,  for  the  creation  of  solid  ice  particles 
directly  from  supersaturated  water  vapor  by  homogeneous  nucleation.  These  functions  are  not 
presented  herein,  since  it  is  the  J(T)  function  for  nucleation  of  ice  from  liquid  water  (i.e.  J^CO) 
that  is  of  interest  herein.  This  function  was  reported  in  a  following  paper. 

This  important  paper  presents  experimental  data  on  nucleation  rates  of  droplets 
condensing  from  supersaturated  cloud  chamber  environments,  and  also,  confirms  that  small 
submicron  and  micron  sized  liquid  droplets  spontaneously  freeze,  in  about  0.01  seconds,  when 
suddenly  formed  and  environmentally  exposed  to  temperatures  at  or  below  -41  °C.  Other 
aspects,  results,  and  observations  of  this  paper  will  not  be  reviewed  herein. 

2.11  Hagen.  Anderson,  and  Kassner.  Jr.  (1981) 

This  paper  [34]  is  a  continuation  of  the  work  presented  in  Anderson,  et  al.,  (1980), 
described  in  2.10.  In  their  1981  paper,  Hagen,  et  al.,  continue  to  analyze  experimental  data  on 
ice  nucleation  to  gain  a  better  understanding  of  ice  nucleation  rates,  as  well  as,  better  estimates  of 
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the  free  energy  of  formation  of  ice  germ  embryos,  under  conditions  of  homogeneous  nucleation, 
i.eM  the  nucleation  of  ice  from  pure  supercooled  water. 


The  authors  discuss  the  important  role  of  homogeneous  nucleation  in  the  ultimate 
objective  of  understanding  heterogeneous  nucleation.  The  authors  present  the  J(T)  function  for 
the  nucleation  rate  of  ice  germs,  or  ice  precursors,  or  ice  embryos,  from  pure  supercooled  liquid 
as  JJT),  in  the  form 


J  is  (T)  =  nv 
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The  parameters  of  this  equation  are  defined  as: 

n'  —  the  number  of  molecules  of  water  in  contact  with  a  unit  area  of  the  ice  germ 


surface, 

v  =  the  volume  of  a  water  molecule  in  ice, 

<ju  s  the  interfacial  free  surface  energy  of  an  ice  and  water  interface, 

nL  =  the  number  of  liquid  water  molecules  per  unit  volume, 

Ag  =  the  activation  energy  for  the  transfer  of  a  water  molecule  across  the  water-ice 
boundary, 

AGls  =  the  increase  in  free  energy  of  the  system  (of  molecules)  of  a  critical-sized 
embryo,  i.e.  an  embryo  that  is  metastable  enough  to  become  an  ice  germ  or  ice 
precursor. 

The  parameters  k  and  h  are  the  Boltzmann  constant  and  the  Planck  constant,  respectively.  The 
expression  for  AG 2s  was  given  as 
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(58) 

where 

ns  =  the  number  of  ice  molecules  per  unit  volume, 

p  =  the  saturated  vapor  pressure  of  water  over  a  plane  surface  of  liquid  water, 
p  =  the  saturated  vapor  pressure  of  water  over  a  plane  surface  of  ice 
This  form  for  JJT)  was  ascribed  to  the  theoretical  developments  of  Turnbull  and  Fisher  [7], 
Dufour  and  Defay  [47],  and  Hobbs  [48]. 

The  authors  note  that  small,  1-20  ^tm  diameter  drops,  were  found  to  freeze  in  time  scales 
of  order  0.01  seconds,  when  exposed  to  ambient  temperatures  at  or  below  about  -40°C. 

By  applying  nucleation  theory  to  fit  the  freezing  populations  in  their  cloud/expansion 
chamber  experiments,  Hagen,  et  al.,  determine  the  “energy  barrier  for  freezing  as  a  set  of 
empirical  fits  for  A g  +  AG^ ,  based  on  the  minimum  temperature  achieved  in  each  experiment. 

Then,  an  averaged  fit  for  A g  +  A was  obtained  by  least  squares,  over  all  the  experiments,  as 

Ag  +  AG'u  =-1.739234 xlO'8 +8.1 157 xlO"2Ir  J  (59) 

in  units  of  joules,  for  temperature  in  Kelvins.  The  authors  present  a  figure  (Fig.  6,  [34])  wherein 
Ag  has  been  separately  estimated  as  a  function  of  T,  given  by 

Ag  =-71.8  x  1  O'20  +  0.3400  xlO-20  T  J  (60) 

where  Ag  has  units  of  joules,  and  T  is  in  Kelvin. 

Note  the  Ag  is  the  same  as  the  parameter  U  in  the  equations  for  J(T)  or  I(T)  that  Mason 
and  co-workers  have  published,  or  the  parameter  A  of  McDonald  [18].  Values  for  Ag,  U  or  A 
range  from  3.3  x  10'20  J  (Mason  [14]),  to  about  9.4  x  10'20  J  (Hagen,  et  al.,  [34]).  Note  that 


25 


McDonald  also  discussed  a  range  of  values  for  Ag  (or,  A)  that  he  estimated  [18]  as  being  from 
about  2.35  x  10  M  to  about  5.5  x  10'2a  J.  Hagen,  et  al.,  discussed  the  differences  in  the  various 
values  of  Ag  and  the  differences  of  interpretation  of  Ag.  The  authors  note  that  Ag,  determined 
from  experiments  in  homogeneous  nucleation,  can  be  used  without  modification  for  application 
to  predicting  heterogeneous  nucleation  events. 

The  authors  incorporate  their  empirical  fit  of  Ag  +  A into  the  expression  for  JJ.T), 

the  homogeneous  nucleation  freezing  nuclei  rate  equation,  to  get 

fl.260.rl05  coo  .1 

!7U(T)  =  7(7)  =  5.92  xlO32  7expj - - - 588.D  (61) 

which  has  units  of  nuclei  per  cubic  centimeter  per  second.  The  authors  show  a  plot  [their  Fig.  7] 
of  log10  JJT)  versus  ambient  temperature,  covering  the  range  of  their  experiments,  and  even 
extrapolating  the  curve  to  greater  ambient  temperatures.  Also  plotted  was  a  J(T)  function 
attributed  to  McDonald  [18],  which  does  fall  within  part  of  their  data  in  the  lower  temperature 
range.  However,  it  is  not  clear  to  the  present  author  that  Hagen,  et  al.,  have  interpreted  the 
experimental  data  appropriately.  The  interpretation,  and  the  evaluation  of  the  J(T)  function  in 
the  present  authors  opinion,  should  be  based  on  the  product  J(T)Vt,  or  its  integral,  which  reflects 
not  only  the  thermodynamic  requirements  of  nucleation,  but  also  reflects  the  size  of  the  droplet 
through  its  volume,  as  well  as  the  time  spent  before  nucleation,  t.  Furthermore,  in  the  present 
author’s  opinion,  when  fit  to  experimental  data,  the  product,  J(T)Vt,  must  reflect  the  active  or 
effective  nucleation  sites,  i.e.,  those  that  lead  to  actual  droplet  freezing,  not  to  just  theoretical 
population  rates.  That  is,  the  nucleation  population-rate  density  function  when  correlated  to 
experiment,  must  subsequently  predict  actual,  effective  freezing  site  population  rates,  not 
theoretical  population  rates. 
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In  summary  of  this  paper,  it’s  importance  includes  the  fact  that  the  authors  analyze  the 
energy  requirements  for  homogeneous  nucleation  freezing  of  small  (generally  less  than  1 
micron)  water  droplets  formed  by  condensation  from  a  supersaturated  vapor  cloud  produced  in  a 
cloud/expansion  chamber.  Their  analysis  shows  that  homogeneous  as  well  as  heterogeneous 
nucleation  freezing  of  water  droplets  can  be  fit  by  empirical  functions  for  Ag  (or  U  or  A)  and 
A Gls  (or  Wc).  Thus,  the  present  author’s  opinion  is  that  when  such  correlations  fit  the  product 
J(T)Vt  to  actual  freezing  data,  the  fit  must  reflect  the  production  rate  of  effective  or  active 
freezing  nuclei,  not  just  the  total  population  rate  of  production. 

2.12  Jensen.  Toon,  and  Hamill  Cl 991) 

This  paper  [35]  represents  another  instance  wherein  homogeneous  nucleation  freezing 
theory  is  modified  to  predict  the  reduction  in  the  freezing  temperature  of  atmospheric  icing 
clouds  by  dissolved  sulfuric  and  nitric  acids  and/or  their  hydrates  in  the  supercooled  droplets. 
The  work  also  relates  directly  to  understanding  polar  atmospheric  cloud  formation. 

The  authors  review  classic  homogeneous  nucleation  theory,  presenting  a  form  of  J(T) 
derived  in  Pruppacher  and  Klett’s  classic  text  Microphvsics  of  Clouds  and  Precipitation  [49], 
namely 


J(T)  =  2NC 


(62) 


in  units  of  nuclei  per  cubic  centimeter  per  second.  The  variables  are: 

Nc  =  number  of  water  molecules  in  contact  with  a  unit  surface  area  of  the  ice  germ, 
embryo,  or  ice  precursor  molecular  cluster. 
pw  =  density  of  liquid  water. 
p(.  =  density  of  (normal)  ice. 

<jiw  =  the  surface  energy  of  the  ice-water  interface. 
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AF*  =  the  energy  that  must  be  overcome  for  the  free  water  molecules  to  become  bound 
to  the  ice  crystal  precursor; 

AF *  is  equal  to  the  difference  between  the  equilibrium  energy  of  the  water  molecule  in  the  liquid 
state,  and  the  energy  of  the  water  molecule  in  the  ice  phase.  This  quantity  has  been  called  the 
phase  change  activation  energy  or  just  the  activation  energy,  and  has  been  approximated  by  the 
energy  of  activation  for  the  displacement  of  water  molecules  in  bulk  water.  AF*  has  been  given 
the  symbol  Ag  by  Hagen,  et  al..  [34],  U  by  Mason  [14],  and  A  by  McDonald  [18]. 

AFg  =  the  work  against  surface  forces  required  to  form  a  critical  ice  germ  of  a  meta¬ 
stable  size,  that  can  become  an  ice  crystal  or  ice  germ,  or  embryo. 

Jensen,  et  al.,  give  AFg  as 

AFg=^icama]  (63) 

where 

ag  =  the  radius  of  the  critical  nucleus  or  ice  germ,  (the  same  as  r ,  before). 

The  expression  for  ag  given  was 


where 

Ma 

K 

L 

T. 


as  =■ 


o i(0 


LmPi  In 


(64) 


Te 

e  J 


=  the  molecular  weight  of  water 
=  the  latent  heat  of  fusion  of  water  (treated  as  positive) 
=  the  melting  temperature  of  ice  (273.16  K) 

=  the  ambient  (supercooled)  temperature 
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These  two  expressions  for  AFg  and  ag  are  the  same  as  that  provided  by,  for  example,  Mason  [14] 


( 1 960),  for  Wc  and  rc,  respectively. 

Jensen,  et  al.,  report  that  ag  is  modified  by  acids  or  acid  hydrates  present  or  in  solution 
with  the  water,  so  that,  based  on  the  theory  of  Pruppacher  and  Klett  [49],  for  this  case 


where 


aUi  =  the  “ice-solute  surface  energy,” 

aa  =  the  activity  of  water  in  the  solution  of  water  and  acid. 

Because  aw  was  reported,  experimentally,  over  a  fairly  wide  temperature  range,  Jensen,  et  al., 
could  determine  the  variation  of  J(T)  over  the  range  of  polar  cloud  temperatures  for  typical 
droplet  sizes.  The  spontaneous  freezing  temperature  of  droplets  now  varies  with  the 
concentration  of  acid  (such  as  H.SOJ  present  in  solution.  Freezing  temperatures  for  solution 

droplets  are  predicted  as  low  as  195K  or  -78C. 

While  this  paper  has  other  important  results  and  conclusions  concerning  the  applicability 
of  modified  classical  homogeneous  nucleation  theory  to  predict  the  freezing  of  the  atmospheric 
aerosol,  these  will  not  be  reviewed  here.  The  main  point  of  the  present  review  is  merely  to 
document  the  form  for  J(T)  used  by  the  authors  and  the  method  used  to  modify  the  function  to 
account  for  the  effects  of  dissolved  material  in  water  to  modify  the  spontaneous  freezing 


temperature  of  very  small,  —1  micron,  aerosol  droplets. 

2.13  Stovanova.  Kashchiev,  and  Kunenova  (1994) 

This  paper  [36]  develops  and  tests  a  method  for  analyzing  and  predicting  the  kinetics  of 
supercooled  droplet  freezing  that  encompasses  both  homogeneous  and  heterogeneous  (seeded) 
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nucleation  freezing.  Classical  homogeneous  nucleation  theory  is  the  starting  point  for  the  kinetic 
freezing  model.  The  aim  of  the  paper  is  “to  propose  a  method  for  experimental  determination  of 
the  nucleation  rate  in  freezing  droplets,  to  employ  the  method  for  obtaining  the  ice  nucleation 
rate,  to  characterize  quantitatively  the  nucleation  activity  of  the  aerosols,  and  to  verify  the 
theoretically  expected  linear  dependence  of  the  nucleation  rate  on  the  aerosol  concentration. 

The  theoretical  considerations  begin  with  the  assumption  that  for  each  resulting  frozen 
drop,  one  ice  germ  or  embryo  was  formed  that  lead  to  the  frozen  drop.  Under  this  hypothesis,  if 
the  population  of  a  sample  of  initially  Na  identical  drops  is  studied  while  being  uniformly 
supercooled  at  a  rate,  q,  then  the  time  rate  of  increase  in  frozen  droplets  should  be  functionally 
related  to  J(T)  (for  both  homogeneous  as  well  as  heterogeneous  freezing  processes).  Stoyanova, 
et  al.,  show  that  time  can  be  eliminated  as  a  variable  of  the  problem  as  follows.  If  dN(t)/dt  is  the 
time-rate-of-change  of  the  number  of  frozen  droplets,  then,  from  the  definition  of  J(t), 

^l^v[N0-N(t)]J(t)  (66) 

dt 

where  v  is  the  volume  of  each  identical  droplet.  For  a  constant  cooling  rate  q,  by  definition 


(67) 


Next,  the  time-rate-of-change  of  the  number  of  drops  freezing  can  be  written 

=  v[Na  -  N(t)]l(t)  (68) 

dt  dT  dt  1  ° 

or,  substituting  for  dT/df  and  rearranging 

_  -~[N0  - N(T)]J(T) .  (69) 

dT  q 

Integrating  this  expression  from  the  initial  temperature,  T0  —  273. 16K,  to  some  supercooled 
temperature,  T,  gives 
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N(T)  =  N0 


1-exp 


--  f  J(T*)dT* 
q  t 


(70) 


Note  that  N(TJ  =  0,  i.e.  there  were  no  frozen  droplets  at  Ta  =  273.16  K. 

Stoyanova,  et  al.,  then  solve  this  equation  for  J(T)  by  first  rearranging  the  terms,  and  then 
differentiating  with  respect  to  T  to  get 


v  dT 


N(T) 


(71) 


This  equation  “shows  that  the  experimental  determination  of  the  temperature  dependence  of  the 
nucleation  rate  J  at  known  droplet  volume  and  constant  cooling  rate  q  reduces  to  finding  the 
temperature  derivative  of  the  experimentally  obtainable  quantity  ln(l  -  N(T)/NJ.”  Also,  this 
procedure  applies  whether  the  droplets  contain  motes  or  foreign  particles  or  not.  Thus,  J(T)  for 
heterogeneous  nucleation/droplet  freezing  events  can  also  be  experimentally  determined.  The 
authors  also  remark  that  this  method  of  determining  J(T)  can  also  be  applied  to  non-steady-state 
nucleation  processes. 

The  authors  next  review  classical  nucleation  theory  and  present  J(T)  in  the  form 

/(r)=A(r)exp{-w*/*r}  (72) 


where 

! 

W*  s  “the  nucleation  work,” 

A(T)  —  “a  kinetic  factor  whose  temperature  dependence  is  usually  weaker  than  that  of 
the  exponential  term.” 

k  is  the  Boltzmann  constant  as  defined  in  a  previous  subsection.  The  authors  note  that  J(T)  is 
really  J(T(t))  or  J(t),  since  the  phase  temperature  usually  depends  on  time.  Application  of  the 
classical  theory  is  valid  for  relatively  slow  temperature  transients  (“...sufficiently  small  cooling 
rate,  q"...),  so  that  J(T)  represents  an  approximation  valid  for  slow  cooling  rates.  That  is,  in 
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other  words,  J(T)  in  the  form  due  to  Volmer  [3],  and  others,  was  derived  for  steady-state 
isothermal  processes.  This  does  raise  questions  about  whether  this  theory  can  be  used  directly  to 
predict  supercooled  droplet  freezing  in  rapidly  accelerating  and  cooling  ducted  spray  flow  fields. 
Clearly,  future  research  on  the  freezing  of  very  rapidly  cooled  water  droplets  is  needed. 
Stoyanova,  et  al.,  present  A(T)  as 

A(T)=ZfNa  (73) 

where 

Z  s  the  Zeldovich  factor 

/  s  the  frequency  of  attachment  of  water  molecules  to  the  molecular  cluster  or  ice 
germ 

Na  s  “the  concentration  of  active  centers  on  which  nuclei  can  be  formed” 

The  Zeldovich  factor,  Z,  was  then  given  as 

Z  =  ty*/3nkTn*2Y2  (74) 

where 

n*  =  “the  number  of  molecules  in  the  nucleus” 

Stoyanova,  et  al.,  give/as 

f  =  n*  fikT /rjvm  (75) 

where 

n*  h  “the  number  of  attachment  sites  of  molecules  on  the  nucleus  surface,” 

r\  =  “the  viscosity  of  the  liquid  around  the  nucleus,” 

vm  =  “the  molecular  volume,”  i.e.  the  volume  of  a  single  molecule, 
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ft  =  is  a  factor  that  accounts  for  the  change  in  liquid  viscosity  near  the  surface  of  the 
molecular  cluster  ((3  <  1). 

The  expression  for  A(T)  that  then  results  is 


Zn]pkT 

Wm 


Stoyanova,  et  al.,  report  that  classical  nucleation  theory  gives 


W*  =  a<jlfvl/Afj. 


and 


(76) 


(77) 


n*  =  2W  */A/t  (78) 

where 

a  =  “a  numerical  shape  factor  for  the  molecular  cluster,  e.g.,  a  =  16tz/3  for  spherical 
nuclei,” 

crf  =  the  specific  (effective)  surface  free  energy  of  the  liquid/nucleus  interface 
(equivalent  to  oSL  or  used  in  other  J(T)  or  I(T)  expressions  presented  herein). 

A/j  =  the  difference  between  the  chemical  potentials  of  the  molecules  in  the  liquid  and 
in  the  ice  germ  or  molecular  cluster. 

The  authors  note  that,  for  heterogeneous  nucleation,  the  value  of  otf  is  less  than  that  for 

homogeneous  nucleation.  So  that,  in  general,  one  can  write 

.  =<E>cr3  (79) 

where 

<T  =  the  specific  surface  free  energy  for  homogeneous  nucleation  (same  as  cru). 

=  0  <  O  <  1;  a  parameter  used  for  quantitative  characterization  of  the  nucleation 
activity  of  aerosol  particles  or  “active  centers.” 
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Stoyanova,  et  al.,  note  that  “different  theoretical  models  give  different  expressions  for  O.”  Thus, 
theoretical  models  such  as  that  discussed  by,  for  example,  Fletcher  [24]  can  be  used  to  estimate 
O  for  the  application  of  J(T)  to  predict  heterogeneous  nucleation  freezing  of  supercooled  water 
drops.  Stoyanova,  et  al.,  report  a  functional  form  for  O  in  the  case  where  a  hemispherical  cap 
shaped  nucleus,  or  ice  germ,  begins  to  form  on  a  flat  surface: 

g,(9)=(2  +  cosgXl-cos8)i  (80) 

4 

This  expression  is  attributed  to  Volmer  [50].  In  this  expression  for  O,  0is  the  angle  of  “wetting” 
of  the  surface  by  the  nucleus.  Values  of  0(0)  would  be 


0(180)  =  1 


for  either  non-wetting  or  homogeneous  nucleation; 

0(90)  =  0.5 
for  “half-wetting”  and; 

0(0)  =  0 


at  “full  wetting.”  Full  wetting  would  seem  to  imply  particle  freezing  at  the  bulk  water  freezing 
temperature,  T  =  273. 16K. 

Stoyanova,  et  al.,  also  provide  thermodynamic  expressions  for  the  evaluation  of  Afi.  For 
the  freezing  of  water,  they  provide 


,  AS  TAT 


(81) 


where 

ASm  =  the  entropy  change  of  melting  (see  equation  (21),  for  example), 

AT  =  Tm  -  T,  the  degree  of  supercooling  experienced  by  the  water  sample, 
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Tm  =  the  melting  temperature  of  ice,  Tm  =  273. 16K. 

By  combining  the  various  expressions  given  so  far,  Stoyanova,  et  al.,  write  J(T)  as 


j(T)  =  ZnsP 


a  eXP' 


A Si  kT3AT 2 


where,  with  a]f  =  <b<73,  J(T)  will  apply  to  both  heterogeneous  and  homogeneous  nucleation 

freezing  events.  For  values  of  some  of  the  parameters  in  the  equation  for  J(T),  the  authors 
provide: 

a  =  1671/3 

<7  =  0.02  J/m2  (or  20  erg/cm2) 

ASm  =  2.65  k  ( k  is  Boltzmann’s  constant) 

v  =  3  x  10'23  cm3 

m 

An  expression  for  the  viscosity  of  liquid  water  was  reported  as 

/  T  y164 

77(r)  =  0.139 - 1  poise  (83) 

\  225  ) 


but,  in  the  temperature  range  of  their  experiments,  the  authors  used 

r\  =  0.005  poise  (constant  rj). 

The  authors  note  that  analysis  of  the  heterogeneous  water  sample  freezing  data  that  they 
had  obtained  then  required  a  value  of  /3  of  about  10*.  Recall  that,  for  Stoyanova,  et  al.,  p  is  the 
factor  that  accounts  for  “the  change  in  liquid  viscosity  very  near  to  the  nucleus  or  ice  germ 
surface.”  The  small  value  for  p  may  “actually  reflect  the  commonly  observed  failure  of  the 
classical  nucleation  theory  to  give  absolute  magnitudes  correctly.”  These  aspects  of  the 
formulation  of  J(T),  due  to  Volmer  [50],  Walton  [51],  and  others,  as  presented  by  Stoyanova,  et 
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al.,  [36],  as  well  as  the  extension  of  the  theory  to  higher  cooling  rates  of  particles,  need  further 
review  and  clarification. 

An  application  of  the  nucleation  theory  was  made  by  Stoyanova,  et  al.,  to  both  seeded 
and  unseed  supercooled  water  freezing.  The  purpose  was  to  evaluate  key  parameters  of  the 
theory,  such  as  aef  O,  6,  Afi,  and  W*  for  both  “seeded”  and  “unseeded”  water  samples.  The 
“unseeded”  water  was  not  pure  enough  to  reflect  homogeneous  nucleation,  being  referred  to  as 
“distilled  water.”  Stoyanova,  et  al.,  provided  tabulated  values  for  the  various  parameters  listed 
above  that  resulted  from  analysis  of  supercooled  water  sample  freezing  due  to  different  types  and 
amounts  or  concentrations  of  “atmospheric  aerosol  particles.”  As  expected,  oif,  O,  6,  Afi,  and 
W*  all  decrease  as  the  concentrations  of  seed  particles  in  the  water  samples  increased.  There 
were  also  changes  in  these  variables  for  different  types  of  seed  particles  or  “active  centers”  as 
they  were  denoted. 

In  summary,  Stoyanova,  et  al.,  like  Fletcher  [24],  describe  a  general  approach  for  the 
application  of  homogeneous  nucleation  theory  to  heterogeneous  nucleation  freezing  processes.  It 
is  clear  that  atmospheric  samples  of  water  collected  from  flight  through  icing  conditions  will  be 
needed  in  the  future  to  begin  the  process  of  modeling  freezing  of  atmospheric  supercooled  water, 
as  well  as  for  the  simulation  of  icing  tests  in  ground  test  facilities.  The  particle  concentrations  in 
these  samples,  as  well  as  classification  of  particulate  types  will  be  needed  for  future  computer- 
based  modeling  studies  of  the  freezing  of  supercooled  water.  Water  droplets  that  impinge  on 
aircraft  surfaces  can  form  thin  run  off  sheets  whose  small  thickness  may  permit  the  use  of 
modified  homogeneous  nucleation  theory,  or  heterogeneous  nucleation  theory,  to  identify  where 
in  the  run  back  process  freezing  begins.  Other  effects  on  the  heterogeneous  nucleation  freezing 
process  in  supercooled  water,  such  as  the  effect  of  different  types  of  seed/mote  particles 
simultaneously  present  in  the  water  samples,  as  well  as  their  concentrations,  were  also  discussed 
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Stoyanova,  et  al.,  which  makes  their  paper  an  excellent  starting  point  for  future  researchers  to 

review  the  theory. 

2.14  Summary  of  the  Overview 

This  section  of  the  report  was  intended  to  provide  a  brief  overview  of  some  aspects  of 
homogeneous  nucleation  theory  applied  to  predict  the  nucleation  or  crystallization  of 
supercooled  liquid  water  droplets.  This  overview  was  not  intended  to  be  an  inclusive  survey  of 
nucleation  theory  and  it  is  not.  It  was  also  the  purpose  of  this  section  to  indicate  how 
homogeneous  nucleation  theory  has  been  extended  to  apply  to  heterogeneous  nucleation  events, 
that  is,  to  instances  where  supercooled  liquid  water  droplets,  of  ordinary  purity,  freeze  due  to 
contact  with  solid  surfaces  or  to  the  presence  of  entrained  foreign  particles  (“motes”)  in  the 
water. 

Because  heterogeneous  nucleation  is  the  predominate  mechanism  by  which  supercooled 
water  drops  freeze  in  natural  and  artificial  icing  environments,  the  use  of  modified  homogeneous 
nucleation  theory  to  predict  droplet  freezing  represents  a  physics-based  freezing  theory, 
compared  to  probability-based  freezing  theories  developed  earlier  by  Levine  [11],  Bigg  [16,17], 
and  applied  by  the  present  author  to  predict  spray  cloud  freezing  in  ground  test  simulations  of 
icing  conditions  [1,2].  The  first-generation  freezing  model  was  based  on  the  Levine-Bigg 
probability  arguments  for  activities  of  the  “freezing  nuclei”  in  mote-induced  freezing.  The 
second-generation  freezing  model  is  based  on  modified  homogeneous  nucleation  theory 
(MHNT)  which  can,  in  principle,  account  for  the  effects  of  mote-types  on  water  particle  freezing 
by  accounting  for  the  “wetting”  characteristics  of  the  mote-types,  as  well  as  for  chemical 
additives  that  change  either  vapor  pressures  or  interfacial  surface  free  energies,  and  for  surface 
geometries. 
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In  the  following  sections,  the  development  of  a  water  particle  crystallization  model  for 
heterogeneous  freezing  is  described.  The  incorporation  of  this  model  into  a  numerical,  one¬ 
dimensional,  multiphase  flow  code  is  discussed  and  results  are  presented  which  were  obtained 
from  predictions  made  of  water  particle  freezing  in  ducted  flows  using  this  code.  Comparisons 
of  these  predictions  are  made  to  the  same  cases  and  similar  results  obtained  with  the  first- 
generation  water  particle  freezing  model  (the  Levine-Bigg  Model)  [1,2]. 

Conclusions  are  drawn  from  the  comparison  of  results  and  recommendations  are  made 
for  further  research  and  development  of  water  droplet  crystallization  models  for  multi¬ 
dimensional  numerical  flow  codes.  The  other  implications  of  modified  homogeneous  nucleation 
theory  are  discussed,  which  include  other  aspects  of  adverse  weather  phenomena,  as  well  as  for 
weather  simulation  in  ground  test  facilities. 

3.0  BEHAVIOR  OF  THE  FREEZING  NUCLEI  RATE  OF  CREATION  FUNCTION,  J(T) 

3.1  J(T)  Utilized  in  Present  Study 

The  first  step  taken  by  the  present  author  to  understand  the  behavior  of  the  function  J(T) 
in  predicting  the  onset  of  freezing  of  submillimeter,  pure,  supercooled  water  droplets  was  to 
reproduce  the  freezing  temperature  versus  droplet  diameter  curve  shown  in  Figure  54,  page  496, 
of  the  paper  of  Langham  and  Mason  [23].  This  curve  was  based  on  their  equation  (4)  of  their 
paper  which  is  provided  as  equation  (43)  of  the  present  report  (which  corrects  the  typographical 
error  in  their  equation  (4)).  The  author  presented  the  freezing  temperature  curve  thus  obtained, 
in  Figure  1,  of  a  previous  report  [1],  as  described  in  [1]. 

In  the  present  study,  the  author  found  it  convenient  to  use  a  simplified  version  of  equation 
(43),  which  is  the  form  first  provided  by  Mason  [14],  shown  in  the  present  report  as  equation 
(36),  page  11.  The  reason  for  using  the  simpler  format  for  the  equation  is  that  by  numerical 
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experimentation,  it  is  obvious  that  the  parameters  having  the  greatest  effect  on  the  value  of  J(T) 
are  uSL,  Lf  and  T.  For  given  values  of  Lf  and  T,  the  function  J(T)  is  very  sensitive  to  values  of  aSL, 
because  its  cube  enters  the  exponential  term  of  J(T).  Therefore,  adopting  a  constant  value  for  Lp 
the  heat  of  melting  of  ice  as 

Lf  =  3.33  xlO5  J/kg, 

and  a  shape  factor  co  evaluated  at 

co  =  23 


the  form  of  the  J(T)  equation  used  in  the  present  study  is  given  by 


J(T)  =  6.9708  xl038r  exp 


2390.2 

2.367  x  107 

T 

T 

_3 

®SL 


(  273.16  Y 


273.16-7 


(84) 


3.2  Properties  and  Behavior  of  J(T) 

The  form  of  J(T)  given  in  equation  (84)  explicitly  includes  cSL  as  a  free  parameter 
meaning  that  crSL  must  be  determined  for  both  homogeneous  as  well  as  heterogeneous  (mote- 
induced)  freezing  processes  which  are  to  be  modeled  in  the  present  study.  In  the  form  given  by 
equation  (84),  assuming,  as  did  Langham  and  Mason  [23],  that  1  micron  diameter  pure  water 
drops  freeze  in  0.6  seconds  at  -41°C,  the  value  for  gsl  was  found  to  be  required  at  . 

o‘SL  =0.0212  J/M2  (85) 

Figure  1  shows  a  plot  based  on  J(T)  of  the  calculated  number  of  critical-sized  active  freezing 
nuclei  for  a  one-micron  diameter  droplet  of  pure  water  as  a  function  of  the  droplet  temperature 
(in  degrees  centigrade).  It  was  assumed,  in  the  present  study,  that  as  a  minimum,  at  least  one 
freezing  nuclei  had  to  be  present  in  the  drop  to  initiate  freezing  at  -41°C.  The  number  of  critical¬ 
sized  nuclei  present  was  evaluated  from  equation  (51)  and  (52)  of  the  present  study,  using  J(T) 
given  by  equation  (84),  under  the  assumption  that  each  value  of  Nc  obtained  was  obtained  for  a 
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fixed  particle  size  (1  ji)  and  at  each  fixed  temperature,  T.  Thus,  under  these  conditions,  equation 
(52)  reduced  to 

NC=J(T)*-D3  *t  (86) 

6 

where  D  =  1  micron,  and  t  =  0.6  seconds.  Note,  that,  the  actual  production  rate  term  for  nuclei 
generation,  that  is,  J(T),  has  a  distribution  with  temperature  as  shown  in  Figure  2.  Thus,  the  use 
of  the  thermodynamic  theory  of  spontaneous  nucleation  to  predict  particle  freezing  cannot  be 
simply  based  on  values  taken  by  J(T),  but  rather,  must  be  obtained  with  J(T)  using  some  criteria 
such  as  that  expressed  by  equations  (51)  and  (52)  in  the  present  study.  As  a  second  point  of 
interest,  Figure  3  is  provided  to  show  the  freezing  nuclei  production  rate  term  J(T)  normalized  by 
the  number  of  molecules  per  cubic  meter  present  in  the  water  sample,  N.  Thus,  Figure  3  shows  a 
plot  of  J(T)/N  versus  the  temperature  of  the  supercooled  water.  As  this  figure  shows,  the 
formation  of  molecular  clusters  containing  between  20-300  molecules,  that  form  the  ice  germs  or 
ice  precursors,  is  still  a  fairly  rare  molecular  process  in  the  water  sample,  even  at  -41°C. 

3.3  Review  of  Dorsch  and  Hacker  Data  for  Heterogeneous  Freezing  of  Submillimeter 
Drops 

In  [1],  the  present  author  selected  the  data  provided  by  Dorsch  and  Hacker  [10]  to  obtain 
a  Levine-Bigg  type  heterogeneous  freezing  function  for  use  in  a  one-dimensional  multiphase 
flow  code.  The  Levine-Bigg  freezing  function  fit  to  the  Dorsch  and  Hacker  data  was 

Tf  =  422.2  +  5.592  In  (D) R  (87) 

where  TF  is  assumed  to  be  the  median  freezing  temperature  of  drops  of  diameter,  D,  in  microns. 

In  degrees  centigrade,  this  function  was 

7>  = -38.6 +  3. 11  In  (Z>)C.  (88) 
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The  Dorsch  and  Hacker  data  set  was  also  selected  in  the  present  study  for  consistency 
with  the  previous  study  [1]  and  for  the  reasons  it  was  selected  in  the  first  study,  as  well  as 
because  this  data  set  falls  in  the  range  typical  of  heterogeneous  freezing  experiments.  The  first 
step  taken  in  the  present  study  to  develop  a  more  physics-based  particle  freezing  function  for 
heterogeneous  freezing,  based  on  modified  homogeneous  nucleation  theory,  was  to  review  the 
Dorsch-Hacker  data  set.  Figure  4  shows  a  set  of  data  plotted  that  came  from  Fig.  8  and  part  of 
Table  I  of  the  Dorsch-Hacker  report  [10].  The  figure  shows  both  the  original,  hemispherical 
“drop”  sizes,  as  well  as  the  equivalent  spherical  droplet  sizes,  Dtq ,  where 

Deq=2~'/iD  (89) 

In  the  present  study,  additional  data  points  were  assembled  from  Table  I,  Fig.  8,  Fig.  10,  and  Fig. 
11  of  the  Dorsch-Hacker  report  which  were  then  combined  into  a  single  figure,  Fig.  5,  of  the 
present  report.  Also  on  Fig.  5  is  a  curve  representing  the  data  utilized  in  [1]  to  obtain  the  Levine- 
Bigg  curve  fit,  equations  (87)  and  (88).  Thus,  the  previously  used  representation  of  the  Dorsch- 
Hacker  data  in  [1]  was  deemed  adequate,  for  defining  a  Levine-Bigg  freezing  function.  Figure  6 
shows  the  data  of  Fig.  5  corrected  to  spherical  diameters.  In  the  present  study,  using  all  of  the 
data  of  Fig.  6  to  obtain  a  new  representation  of  a  Levine-Bigg  freezing  function,  a  new  curvefit 
equation  was  obtained  for  particle  heterogeneous  freezing  temperature,  given  by 

Tf  =  -38.3  +  3. 19  In  (D)  in  degrees  C  (90) 

This  function  is  plotted  in  Fig.  7,  and  is  essentially  the  same  as  that  obtained  previously  [1],  i.e., 
it  compares  favorably  to  equation  (88). 

Therefore,  either  of  the  Levine-Bigg  type  representations  of  the  Dorsch-Hacker  data  are 
assumed  accurate  enough  for  expressing  freezing  behavior  of  the  Dorsch-Hacker  experiments  on 
heterogeneous  particle  freezing. 
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The  Levine-Bigg  type  of  freezing  function  for  heterogeneous  freezing  of  water  particles 
is  considered  to  be  a  first-order  or  lowest-level  freezing  model  from  the  standpoint  of  the  physics 
of  the  process. 

3.4  Freezing  Function  Based  on  Modified  Homogeneous  Nucleation  Theory  (MHNT) 

A  second  order,  or  higher-level  freezing  model,  was  developed  next  in  the  present  study 
based  on  the  following  criterion  and  equations: 

(1)  Criterion  for  initiation  of  freezing  of  a  water  particle; 
with  the  number  of  active  freezing  nuclei  being  given  by 

t 

Nc(T,D,t)=\j{T,D,t)V{D,t)dt  (91) 

the  freezing  criterion  is 

Nc(T,D,tF)=  1.0  (92) 

That  is,  in  the  present  study,  it  is  assumed  that  at  least  one  active  freezing  nuclei  must  be  created 
in  the  drop  of  volume,  V,  in  the  time  period  T,  whereT  =  tF  —  ta ,  to  initiate  freezing.  The  times 

are  defined  herein  as 

to  =  the  time  when  the  droplet  temperature  first  reaches  the  bulk  water  freezing 
temperature,  i.e.  273. 1 6  K 

tF  -  the  time  when  the  drop  begins  to  freeze,  i.e.,  by  definition,  when  Nc=  1.0 

(2)  It  is  assumed  that  the  homogeneous  nucleation  rate  function,  J(T),  can  be  used,  as 
modified,  for  heterogeneous  freezing  predictions.  That  is,  the  freezing  nuclei  rate  equation, 
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J(T,D,t),  is  the  same  function  as  the  homogeneous  nucleation  rate  equations  given  by  equation 
(84),  except,  that  cSL  is  now  a  function  of  drop  diameter;  that  is,  cSL  -  <JSL  (D).  Thus,  J(T,D,t)  for 
heterogeneous  freezing  is  assumed  given  by 


j{T,D,t)=  6.9708  xl038  rexp 


2390.2 

2.367 JtlO7 

T 

T 

(  273.16  Y 
273.16-T 


(93) 


(3)  It  is  assumed  that  the  specific  surface  free  energy  for  heterogeneous  freezing  is  a 
function  of  drop  volume,  hence  a  function  of  the  number  and  type  of  active  freezing  nuclei 
present  in  the  droplet  volume.  For  simplicity,  at  present,  therefore,  it  was  assumed  that  oSL  could 
be  represented  by  a  smooth  function  of  the  droplet  diameter,  D,  or, 

Osl  =  aSL  ( D )  (heterogeneous  freezing) 

From  a  practical  standpoint,  therefore,  the  specific  surface  free  energy,  oSL(D),  for  the 
heterogeneous  nuclei  (water  plus  mote)  responsible  for  the  initiation  of  droplet  freezing  in  the 
Dorsch-Hacker  experiments  can  be  determined  for  representative  drop  sizes  by  matching  the 
criterion  given  by  equation  (92)  for  each  drop  size,  at  its  median  freezing  temperature.  To 
implement  this  process,  and  to  develop  a  smoothly  varying  oSL  function  of  D,  the  median 
freezing  temperature  of  each  size  drop  was  assumed  to  be  given  by  equation  (90).  Thus,  oSL(D) 
was  determined  by  requiring  Nc  (TFP  DP  tFi)  to  be  unity  for  a  selected  set  of  freezing  temperatures 
predicted  by  equation  (90).  Thus,  if  TFi  is  the  median  freezing  temperature  of  droplet  with 
diameter,  Dv  then 

N,(Tri.D„tK)=^0  =  )j(TFl,Dnl)V(D„t)dt  (94) 

lOi 

TFi  and  V(DP  t)  were  assumed  constant  during  an  assumed  freezing  interval  given  by 

t.  =  tFi  -tol  =  1.0  second. 
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Then  NC(DP  TF.)  could  be  approximated  by 

W[(7-fl,Dr)=1.0  =  7(rn,D,)V(D,)(l.0)  (95) 

With 

V(D1.)=-D13  (96) 

6 

the  freezing  criterion  actually  applied  to  determine  GSL(D)  was 

Nc  (T F, ,  A  ) = 1  -0  =  J  (rF, .  Di  )|  A3  (97) 

Figure  8  shows  N/T„  D.)  plotted  as  a  function  of  T  for  drop  sizes  ranging  from  1-10000 
microns.  The  line  Nc  =  1.0  cuts  each  Nc  curve  at  the  median  freezing  temperature  of  that 
respective  droplet  size.  For  example,  for  the  5  micron  droplet,  the  line  Nc  =  1.0  cuts  the  5  micron 
Nc  curve  where  T  =  -33°C.  By  inspection  of  Fig.  7,  it  can  be  seen  that  this  corresponds  to  the 
point  on  the  curve  fit  for  the  heterogeneous,  spontaneous  (median)  freezing  temperature  where  5 
micron  droplets  freeze  at  TF  —  -33°C  (the  ordinate  value  of  TF  for  In  (5)  =  1.609). 

Each  curve  plotted  on  Fig.  8  was  obtained  by  trial  and  error  evaluation  of  <JSL,  changing 
the  value  of  gsl  for  each  drop  size,  until  the  Nc  curve  generated  for  that  drop  size  passed  through 
the  appropriate  freezing  temperature  when  Nc  =  1.0. 

With  the  determination  of  each  successful  value  of  GSL,  for  each  a’ priori  chosen  drop  size, 
the  data  in  Fig.  9  was  constructed.  Shown  in  Fig.  9  are  two  curve  fits  for  aSL  that  were  obtained 
by  the  same  method  to  deduce  the  variation  of  aSL  for  heterogeneous  freezing.  The  first  attempt 
produced  the  distribution  denoted  GSU(D).  A  second  attempt  to  create  a  oSL  distribution  in  the 
present  study  corresponds,  in  Fig.  8,  to  the  curve  osu(D).  The  differences  between  Gsu  and  Gsl2 
result  from  a  later,  more  thorough  reconstruction  and  analysis  of  the  Dorsch-Hacker  data  by  the 
author.  The  forms  of  the  curvefit  equations  for  Gsu  and  Gsu  are 
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asu  (D)  =  0.02045  -  0.00044528  a:  (in  (D))-  7.2527  *  10'5  (ln(Z)))2 


(98) 


and 

oSL1  =  0.020228  - 0.00044755  *  (ln(D))- 7.908 1  x  10"5  x  (in (D)f  (99) 

These  two  representations  for  aSL  were  kept  in  the  present  study  to  investigate  the 
sensitivity  of  the  predicted  freezing  behavior  to  variations  in  the  values  of  aSL,  for  the  same 
water  sample  data  set.  <Jsa  is  the  recommended  curvefit  because  it  is  believed  more  accurate. 

The  differences  in  predicted  submillimeter  water  particle  freezing  behavior  in  ducted 
flow  when  computed  using  the  Levine-Bigg  theory,  (7SU,  and  oSL2,  will  be  shown  in  the  next 
section. 

4.0  ANALYSIS  AND  EVALUATION  OF  THE  WATER  PARTICLE  FREEZING 
MODEL  BASED  ON  MODIFIED  HOMOGENEOUS  NUCLEATION  THEORY 
4.1  Introduction  to  the  Approach 

Because  of  a  lack  of  detailed  data  on  the  freezing  of  convected,  submillimeter, 
supercooled  water  droplets  in  ducted  flows,  the  analysis  and  evaluation  of  the  new  model  for 
heterogeneous  freezing  must  be  based  on  comparing  model  results  with  model  results. 
Therefore,  the  predicted  results  of  a  calculation  of  water  particle  freezing  in  a  ducted  flow,  made 
based  on  the  Levine-Bigg  freezing  model  [1],  will  be  compared  to  the  corresponding  results 
obtained  based  on  the  modified  homogeneous  nucleation  theory  (MHNT)  presented  in  the  last 
section.  The  results  obtained  with  the  Levine-Bigg  model  were  made  with  a  one-dimensional, 
multiphase  flow  code  described  in  references  52,  53  and  modified  and  updated  in  the  study 
reported  in  [1].  The  flow  code  solves  the  fully-coupled  mass,  energy,  and  momentum 
conservation  equations  for  a  dilute,  multiphase  flow  of  water  particles  entrained  in  a  ducted 
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airflow  [52,  53].  The  solution  procedure  is  a  Runge-Kutta  fifth-order  numerical  integration 
scheme  [54,  55]  that  requires  a  set  of  input  duct  flow  initial  conditions  as  well  as  the  duct 
geometry. 

The  code  was  denoted,  in  the  author’s  first  report  on  water  particle  modeling  [1],  as  the 
“AEDC1DMP”  code  (pg.  34,  [1]).  This  code  designation  will  be  retained  in  the  present  study, 
meaning  AEDC  one-dimensional,  multiphase  flow  code.  It  will  be  made  clear  when  the  results 
of  this  code  refer  to  predictions  made  based  on  the  Levine-Bigg  (LB)  type  freezing  model  or  the 
modified  homogenous  nucleation  theory  (MHNT)  model  for  particulate  freezing. 

4.2  Implementation  of  the  Particle  Freezing  Model 
4.2.1  Freezing  Criterion 

In  implementing  the  MHNT  model  in  AEDC1DMP,  the  transformation  of  the  calculation 
of  the  time  integral  for  NJT,  D,  t)  to  an  integral  along  the  flow  with  respect  to  the  axial  flow 
coordinate,  was  done  as  follows.  The  time  integral  was  rewritten  as 

N  (jtD,x)~  —  (100) 

i  m 

where 

dx  =  U{x)dt  (101) 

and  U(x)  is  the  particle  velocity  at  x.  The  lower  limit  on  the  integral,  xo,  is  the  axial  location  in 
the  flow  where  the  given  sized  droplet  first  reaches  (or  cools  down  to)  a  temperature  of  32°F  (0 
C  or  273.16  K).  Thus,  droplet  temperatures  are  monitored  and  when  a  given  size  drop 
temperature  reaches  32°F,  the  evaluation  of  the  integral  is  begun.  However,  the  integral  was 
actually  implemented  in  the  code  by  a  running,  finite  summation 

Nfi,D,Xi)=±  002) 
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where  /,  V,  and  if  are  mean  or  representative  values  of  J,  V,  and  U  over  the  spatial  integration 
step  interval,  <£t  ,  at  the  axial  location,  xk,  given  by 

x  =  xk  =  *,+£  dx,  (103) 

I 

The  integration  step  intervals,  dx.,  are  controlled  by  the  numerical  integration  scheme  to  control 
numerical  (truncation)  errors.  Note,  that  when  x  reaches  a  limiting,  input  preset  value,  x^,  then 
the  AEDC1DMP  code  terminates  the  calculation  (integration)  process  for  the  entire  flow. 

The  implementation  of  the  Levine-Bigg  type  freezing  function  in  AEDC1DMP  was 
discussed  previously  [1,2].  However,  to  review  the  process,  the  implementation  was  as  follows. 
In  the  course  of  computing  the  flow,  the  temperatures  of  the  liquid  drops  are  monitored.  When  a 
drop  of  diameter  D  cools  to  a  temperature  at  or  below  the  corresponding  median  freezing 
temperature  TLf ,  for  that  drop  size,  the  droplet  is  considered  to  have  begun  the  freezing  process, 

that  is,  the  freezing  procedure  is  initiated. 

4.2.2  Model  of  the  Freezing  Process 

Whether  the  initiation  of  freezing  is  triggered  by  a  criterion  given  by  the  Levine-Bigg 
model  or  by  the  modified  homogeneous  nucleation  theory  (MHNT)  model,  once  freezing  has 
been  initiated,  the  computation  of  the  drop  freezing  is  the  same.  The  steps  are 

(1)  the  supercooled  liquid  drop  is  assumed  to  become,  instantaneously,  a  slushball  or 
mixture  particle  containing  both  liquid  water  and  ice  at  the  normal  melting 
temperature  of  ice  (492°R,  32°F,  or  0°C). 

(2)  The  amount  of  the  mixture  particle  that  is  ice  is  determined  from  an  energy  balance 
for  the  particle  based  on  conditions  just  before  freezing  was  initiated  and  just  after 
freezing  was  initiated.  Since  the  mixture  temperature  is  assumed  at  492°R,  the  only 
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free  variable  that  can  be  determined  from  the  energy  balance  is  the  fraction  of  liquid 
that  remains  after  freezing  is  triggered,  given  by 

a,  =l-|M*92-rJ  <104> 

hf 

where 

aF  =mass  fraction  of  drop  or  mixture  particle  that  remains  liquid  at  492°R.  (hence, 

(1  -  aF)  of  the  particle  has  become  ice  at  492°R). 

CL=  specific  heat  of  liquid  water 
Hr  e  heat  of  fusion  of  liquid  water 

7  =  freezing  temperature  of  the  supercooled  liquid  water  droplet. 

Details  of  the  energy  balance  are  provided  in  [1]. 

(3)  The  density  of  the  mixture  particle  is  now  the  mean  density  of  a  two-phase  system 
wherein  the  phases  cannot  occupy  the  same  volume.  Hence,  the  mean  density  of  the 
mixture  particle  is  given  by 

p-  Ej£l _  (105) 

aFp,+(\-aF)pL 

where 

pL  =  density  of  liquid  water 
p,  =  density  of  ice 

The  specific  heat  of  the  particle  is  also  redefined  as  a  mass-average  specific  heat 
given  by 

Cr=ar  Cn  +(l-aF)CPj  (106) 

where 
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CPl  =  specific  heat  of  liquid  water 
CP/  =  specific  heat  of  ice 

(4)  The  size  of  the  mixture  particle  is  then  given  by  requiring  that  the  total  mass  of  the 
particle  remains  constant  during  the  (assumed  instantaneous)  transition  of  the  particle 
from  supercooled  liquid  drop  to  a  mixture  particle  of  water  and  ice  at  492°R.  Thus, 

PtD;  =  Pl  7  A’  (107) 

0  o 


or 


Db 


(108) 


where 

D  =  the  new  diameter  of  the  mixture  particle 

DLf  s  the  diameter  of  the  supercooled  liquid  drop  just  before  freezing  is  initiated. 

(5)  Thereafter,  the  mixture  “drop”  is  treated  as  any  other  droplet  in  the  flow.  Heat,  mass, 
and  momentum  exchanges  with  the  air  flow  are  computed  in  the  same  way  as  that  for 
the  all  droplets.  However,  as  each  amount  of  heat  is  extracted  from  the  particle,  the 
amount  of  liquid  remaining  in  the  droplet  is  reduced  ( cxF  is  reduced),  and  a  new  mean 
density,  and  particle  diameter  are  re-calculated.  In  the  recalculation,  the  mass  loss  or 
gain  of  the  particle  by  mass  transfer  is  taken  into  account,  as  well  as  energy  lost  or 
gained  by  the  mass  transfer  process.  The  amount  of  new  ice  formed  is  also 
calculated.  The  details  of  the  mass  and  energy  balance  during  freezing  of  the 
evaporating,  two-phase  particle  are  given  in  Appendix  B. 

When  the  entire  particle  has  become  ice,  preaches  zero,  and  the  particle  is  now 
completely  frozen. 
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During  the  freezeout  process  for  the  mixture  particle,  the  particle  temperature  is 
held  constant  at  492°R. 

(6)  Once  the  particle  completely  freezes,  its  density  becomes  that  of  ice,  as  well  as  its 
specific  heat.  From  this  instant  or  location  in  the  duct  flow,  the  particle  undergoes 
mass,  energy,  and  momentum  exchanges  with  the  air  flow  as  an  ice  particle.  It  can 
gain  or  loose  mass  by  sublimation  based  on  the  partial  pressure  of  ice  and  the  air  flow 
specific  humidity. 

5.0  EVALUATION  OF  AEDC1DMP  CODE  ON  REPRESENTATIVE  DUCT  FLOW 
CASES 

5. 1  Baseline  Case  of  Ducted.  Two-Phase  Flow 

The  base  case  computed  for  a  comparison  purpose  was  a  ducted  flow  with  a  “single” 
water  spray  station  at  its  inlet. 

In  the  code  AEDC1DMP,  the  water  is  assumed  to  enter  the  airflow  at  up  to  ten  different 
injection  stations  that  tire  separated,  axially  in  the  duct.  Each  injection  station  is  capable  of 
putting  in  a  given  amount  of  water,  in  a  given  sized  water  particle  with  its  specified  velocity  and 
temperature.  When  all  of  the  injection  stations  are  bunched  close  together,  axially,  in  the  duct, 
they  can  be  used  to  model  a  single  spray  station  with  a  spray  droplet  size  distribution 
characterized  by  ten  discrete  drop  sizes.  This  was  the  approach  taken  in  the  previous  study  [1]  to 
examine  the  freezing  of  two-phase,  ducted  flows  representative  of  both  full-scale  and  research- 
scale  icing  test  facilities.  Note,  that  to  illustrate  and  emphasize  the  ducted-flow  freezing 
processes  of  supercooled,  submillimeter  water  particles,  the  air  flow  inlet  total  temperature 
utilized  was  significantly  lower,  at  460°R,  than  is  typical  of  many  icing  tests,  where  an  airflow 
inlet  total  temperature  of  486°R  is  more  likely.  In  any  case,  for  consistency  with  previously 
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obtained  results,  the  same  airflow  and  water  inlet  conditions  have  been  specified  and  are 
provided  below. 


Table  1.  Duct  Air  Flow  Inlet  and  Water  Spray  Conditions  for  Baseline  Test  Case  1 


Air  Inlet  Conditions 


Total  Pressure,  psfa  1 123.48 

Total  Temperature,  °R  460°R 

Velocity,  ft/s  20 

Mach  Number  0.0189 

Inlet  Relative  Humidity  .percent  15.78 


Water  (Particle)  Input  Conditions 

Injection  Station  No. 

(ft) 

Particle  Diameter 
(microns) 

Particle  Temperature 

(°R) 

Particle  Velocity 
(ft/s) 

Load  Factor 
(FL) 

1.  0.0 

530 

46 

6.85  E-05 

iPPPPBFgiIiT1^|1 

10 

530 

46 

2.06  E-04 

15 

530 

46 

4.80  E-04 

20 

530 

46 

3.62  E-04 

30 

530 

46 

2.06  E-04 

6.  0.05 

40 

530 

46 

3.43  E-05 

50 

530 

46 

1.37  E-05 

60 

530 

46 

1.37  E-06 

9,  0.08 

80 

530 

46 

1.37  E-07 

10,  0.09 

100 

530 

46 

1.4E-08 

Total  FL  =  1.37E-U3 

FL  =  lbm  of  water/lbm  of  dry  air  injected  at  a  given  injection  station 


The  duct  geometry  used  for  this  reference  or  baseline  case  is  shown  in  Fig.  10.  As  is 
clear  from  inspection  of  Fig.  10,  and  as  discussed  in  [1],  this  duct  geometry  is  not  representative 
of  full-scale  icing  test  facilities.  Rather,  it  was  chosen  in  the  earlier  study  [1]  because  it  was 
related  to  a  similar  geometry  of  interest.  For  consistency  in  comparing  results,  therefore,  this 

geometry  was  also  retained  in  the  present  study. 

5  2  Comparison  of  Predicted  Results  of  Baseline  Case  for  Three  Particulate  Freezing 


Models 

Figure  1 1  shows  predicted  temperatures  of  10  micron  water  droplets,  as  they  vary  with 
axial  distance  along  the  flow,  for  three  heterogeneous  freezing  models.  The  models  utilized  were 
the  Levine-Bigg  model,  given  specifically  by  equation  (87),  and  the  modified  homogeneous 
nucleation  theory  (MHNT)  using  osu  and  csu  for  the  specific  surface  free  energy  specifications. 


51 


As  can  be  seen,  the  comparison  is  extremely  good.  The  results  obtained  with  the  cxSL  function  asu 
gives  a  closer  agreement  with  the  Levine-Bigg  model  than  does  asu.  Since  it  is  thought  that  <7sa 
was  obtained  from  a  more  thorough  data  analysis  than  asu,  the  better  agreement  of  the  predicted 
results  with  asu  is  expected,  although  such  good  agreement  in  both  quality  and  quantity 
displayed  in  Fig.  1 1  is  remarkable. 

On  the  basis  of  the  results  shown  in  Fig.  11,  and  found  to  be  consistently  the  same  for  all 
droplet  sizes,  it  was  decided  to  utilize  the  function  osJD)  for  the  heterogeneous  freezing  model 
for  all  further  calculations  made  in  the  present  study. 

Figure  12  shows  a  comparison  of  predicted  droplet  temperatures  made  with  two  versions 
of  AEDC1DMP,  one  version  with  the  Levine-Bigg  (LB)  freezing  model  implemented,  and  the 
other  version  with  the  modified  homogeneous  nucleation  theory  (MHNT)  model  implemented, 
based  on  osn.  The  agreement  between  the  first  (LB)  and  second  order  (MHNT)  models  is  good. 
An  interesting  difference  in  the  predicted  results  occurs  for  the  5  micron  droplets.  The  Levine- 
Bigg  model  results  in  the  5  micron  droplets  freezing  at  the  location  in  the  duct  flow  where  their 
temperature  has  reached  about  427°R.  (The  freezing  occurs  so  rapidly,  it  occurred  between  print 
step  sizes  in  the  program  and  is  not  shown,  therefore,  in  the  plot,  Fig.  12.)  The  drops  have 
shrunken,  by  evaporation,  to  about  2  micron  in  diameter  at  this  location.  On  the  other  hand,  the 
modified  homogeneous  nucleation  theory  model  predicts  that  the  5  micron  drops,  shrunken  to  2 
micron,  do  not  freeze  within  the  domain  of  the  calculation,  because,  even  though  the  drops  have 
reached  a  freezing  temperature,  their  volume  is  too  small  to  generate  at  least  one  active  freezing 
nuclei.  In  other  words,  the  criterion 

Nc=\jVdt  =  1  (109) 
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was  not  satisfied  for  these  small  drops  throughout  the  calculated  length  of  flow.  Thus, 
significant  potential  differences  could  arise  in  predicted  freezing  behaviors  made  with  a  Levine- 
Bigg  model  and  with  a  modified  homogeneous  nucleation  theory  model,  for  certain 
circumstances. 

The  predicted  temperature  distributions  for  all  of  the  other  drop  sizes  appear  similar  or 
comparable,  based  on  the  two  different  freezing  models. 

(The  calculated  temperature  distributions  shown  in  Fig.  12  also  indicate  that  the  particles 
took  longer  (more  distance)  to  freeze  for  the  MHNT  model  predictions,  than  in  the  Levine-Bigg 
model  predictions,  than  in  the  Levine-Bigg  model  predictions.  This  was  not  due  to  any 
differences  caused  by  the  freezing  models,  but  due  to  improvements  made  in  AEDC1DMP  to 
integrate  the  heat  and  mass  transfer  processes  of  the  freezing  particles  in  a  more  recent  version  of 
AEDC1DMP  in  which  the  MHNT  model  was  implemented.) 

Figure  13  shows  a  comparison  of  the  predicted  droplet  sizes  from  AEDC1DMP  with  the 
MHNT  model  implemented,  for  the  ten  different  drop  sizes  as  they  move  along  the  flow.  Drop 
freezing  shows  up  in  these  curves  by  sudden  small  increases  in  drop  size  (a  kink)  in  the  curves. 
Clearly,  the  5  micron  particle  curve  does  not  show  a  sudden  increase,  indicating  that  the  drop 
stays  liquid  along  the  entire  flow. 

5.3  Calculations  of  Flow  in  a  Representative  Icing  Test  Facility 
5.3.1  Introduction 

The  duct  geometry  of  baseline  case  calculated,  with  results  shown  in  Figures  10-13  was 
not  representative  of  full  scale  icing  test  facility  duct  geometries.  Therefore,  a  second  test  case 
was  prepared  and  calculated  wherein  the  duct  or  wind  tunnel  geometry  and  scale  represents 
typical  icing  research  or  icing  test  wind  tunnel  geometry. 

The  results  from  computations  of  test  case  2  are  shown  below. 


5.3.2  Nominal  Test  Conditions  in  a  Full-Scale  Icing  Facility 


The  duct  or  wind  tunnel  geometry  assumed  for  test  case  2  is  shown  in  Figure  14.  This 
wind  tunnel  like  configuration  has  a  large  inlet  and  a  steep  contraction  section,  to  minimize  flow 
disturbances.  Icing  spray  nozzles  are  assumed  to  be  put  in  the  plane  of  the  large  inlet.  The  test 
section  is  located  at  x  =  46  ft,  so  that  the  state  of  the  water  particles  injected  at  the  inlet  of  the 
wind  tunnel  should  be  in  kinetic  and  thermal  equilibrium  with  the  air  flow  by  the  time  they  reach 
the  test  section.  The  tunnel  has  a  large  contraction  ratio  (inlet  flow  area  divided  by  test  section 
flow  area),  hence,  the  range  of  permissible  inlet  air  flow  velocities  is  small,  ranging  from  near  0 
to  about  45  ft  per  second.  At  the  higher  inlet  air  velocities,  the  test  section  has  reached  near- 
choking,  or  sonic  flow  conditions.  Generally,  icing  tests  are  conducted  at  air  flow  or  flight- 
simulating  speeds  in  the  range  of  a  few  hundred  miles  per  hour.  Therefore,  a  nominal  set  of  inlet 
test  conditions  was  defined  for  the  representative  calculation  made  with  AEDC1DMP.  The 
nominal  test  conditions  are  listed  below. 


Table  2.  Nominal  Tunnel  Inlet  Test  Conditions  for 
Icing  Wind  Tunnel-Test  Case  2 

Air  Inlet  Conditions 


Total  Pressure,  psfa 

2074 

Total  Temperature,  °R 

500 

Velocity,  ft/s 

20.8 

Mach  Number 

0.0188 

Relative  Humidity,  % 

54.7 

Water  (Particle)  Input  Conditions 


Injection  Station  No. 
(ft) 

Particle  Dia. 
(microns) 

Particle  Velocity 
(ft/s) 

Load  Factor 
(FL) 

_ (a)  1  (0-0) _ 

■■KMfgi 

634 

15.6 

0.00018 

— . . 

500 

634 

15.6 

0.00018 

Note:  The  case  was  computed  first  (case  2a)  with  a  single  particle  size  of  50  microns  to 
represent  the  spray  cloud,  then  computed  (case  2b)  with  a  single  particle  size  of  500  microns.  In 


each  computation,  the  50  or  500  micron  particle  represented  the  mean  volumetric  diameter  of  the 
spray  cloud  particle  size  distribution  at  the  inlet. 

Some  selected  results  from  computing  these  flow  cases  are  shown  in  Figures  15-17. 
Figure  15  shows  the  gas  velocity  together  with  the  velocities  of  a  50  micron  particle  and  a  500 
micron  particle.  The  calculation  shows  that  the  larger  particle  did  not  reach  velocity  (kinetic) 
equilibrium  with  the  air  flow  at  the  test  section.  On  the  other  hand,  the  50  micron  particle  did 
achieve  a  velocity  equal  to  the  gas  velocity  at  the  test  section  location. 

Figure  16  shows  similar  results  for  the  gas  temperature  and  the  temperatures  of  the  two 
different  sized  particles.  In  the  case  of  the  500  micron  particle,  its  temperature  is  above  the  gas 
temperature  by  about  10  degrees,  R,  at  the  test  section  location  due  to  thermal  lag,  however,  the 
50  micron  particle  temperature  has  fallen  quickly,  through  evaporation  and  cooling,  to  its  wet- 
bulb  temperature  about  5  degrees,  R,  below  the  gas  temperature.  The  effects  of  air  flow  inlet 
humidity  on  droplet  cooling  and  thermal  non-equilibrium  have  been  well-documented  elsewhere, 
for  example  [53],  and  will  not  be  discussed  herein.  However,  it  is  possible  to  increase  the  inlet 
air  flow  humidity  to  increase  the  drop  wet  bulb  temperature  at  the  test  section,  to  be  equal  to  the 
gas  temperature.  Large  droplet  thermal  lag  effects  could  also  be  reduced,  by  reducing  inlet 
humidity.  Thus,  opposite  humidity  modifications  are  required  for  large  and  small  drops  to 
account  for  thermal  non-equilibrium  effects  at  the  test  section. 

Finally,  Figure  17  shows  the  changes  in  the  droplet  diameters  as  the  water  particles  pass 
down  the  duct.  The  changes  in  droplet  diameter  are  greatest  where  the  rate  of  evaporation  is 
greatest,  in  the  duct  inlet  throat  region.  (Note  that  kinks  or  slope  discontinuities  in  the  plotted 
curves  are  due  to  the  data  output  frequency  or  axial  spacing,  not  to  discontinuities  in  computed 
results.) 


Test  Case  2,  therefore,  provided  some  nominal  results  which  show  how  AEDC1DMP  can 
be  used  for  assessment  of  flow  variations  and  flow  quality  typical  of  full-scale  icing  test  or 
research  wind  tunnels.  In  this  test  case,  neither  the  50  micron  drop  nor  the  500  micron  drop  were 
near  their  spontaneous  (median)  droplet  freezing  temperatures.  The  next  test  case  computed  will 
be  one  that  represents  a  full-scale  icing  test  with  supercooled  droplet  freezing. 

5.3.3  Test  Conditions  in  a  Full-Scale  Icing  Facility  with  Supercooled  Droplet 

Freezing 

By  trial,  computation,  and  output  analysis,  a  set  of  wind  tunnel  inlet  conditions  was  found 
that  produced  supercooled  droplets,  slushball  or  mixture  particles,  and  ice  particles  in  the  test 
section  located  of  the  flow.  The  set  of  inlet  conditions  found  is  listed  below  in  Table  3. 

Table  3.  Tunnel  Inlet  Test  Conditions  for  Icing  Wind  Tunnel  Flow  with 
Water  Particle  Freezeout  -  Test  Case  3 

Air  Inlet  Conditions 


Total  Pressure,  psfa 

2076 

Total  Temperature,  °R 

475 

Velocity,  ft/s 

40 

Mach  Number 

0.037 

Relative  Humidity,  % 

26 

Water  (Particle)  Input  Conditions 


Injection  Station  No. 
(ft) 

Particle  Dia. 
(microns) 

■■KhH 

Particle  Velocity 
(ft/s) 

Load  Factor 
(FL) 

1  (0.0) 

5 

495 

16 

9.06  E-06 

10 

495 

16 

2.73  E-05 

3  (0.04) 

15 

495 

16 

6.35  E-05 

20 

495 

16 

4.79  E-05 

5  (0.08) 

30 

495 

16 

2.73  E-05 

6  (0.10) 

495 

16 

4.54  E-06 

7  (0.12) 

50 

495 

16 

1.31  E-06 

■KKiXEnMi 

60 

495 

16 

1.81  E-07 

80 

495 

16 

1.8  E-08 

100 

495 

16 

Total  FL  =  1.82  E-04 

S 


As  the  input  data  in  Table  3  indicate,  an  inlet  spray  cloud  with  a  10  drop  size  distribution 
was  input  at  the  wind  tunnel  entrance  station.  The  inlet  air  temperature  and  the  water  spray 
temperatures  are  rather  low,  the  water  being  only  slightly  above  freezing  and  the  air  at  below 
freezing  temperature  (15°F).  This  combination  of  inlet  conditions  lead  to  predicted  freezeout  of 

some  of  the  drop  sizes  in  the  wind  tunnel  flow  at  the  test  section  station. 

Figures  18  show  some  of  the  calculated  results  obtained  with  AEDC1DMP  utilizing  the 
modified  homogeneous  nucleation  freezing  model  for  predicting  the  freezing  of  supercooled 
water  particles.  Figure  18  shows  the  air  flow  Mach  number  distribution  from  duct  entrance  to 
the  test  section.  Note  that  this  case’s  inlet  conditions  lead  to  rather  high  test  section  Mach 
number,  about  0.69.  This  is  about  twice  the  nominal  or  typical  value  for  icing  tests  in  full-scale 
facilities.  The  higher  value  of  air  flow  Mach  number  ensured  that  air  flow  static  temperatures 
would  be  low  enough  to  induce  particle  freezeout,  with  a  freezing  model  calibrated  to  the 
Dorsch-Hacker  data  set  [10]. 

Figure  19  shows  the  predicted  variations  in  drop  or  particle  sizes  from  entrance  to  test 
section.  Note  that  the  initially  5  micron  particle  has  disappeared,  by  evaporation,  by  about  4.5- 
5.0  feet  from  the  duct  entrance.  This  indicates  that,  for  similar  inlet  conditions,  a  real  water 
spray  cloud  with  similar  drop  size  distributions  would  experience  all  of  the  drops  smaller  than, 
say,  10  microns  evaporating  quickly  in  the  inlet.  The  smallest  surviving  droplet  at  the  test 
section  would  probably  be  in  the  range  from  6-8  microns.  These  results,  of  course,  depend  on 
the  inlet  conditions,  especially  water  injection  temperature  and  inlet  air  flow  humidity. 

Figure  19  also  shows  size-jumps  in  the  15  to  80  micron  particles  which  indicates  that 
these  water  particles  have  begun  to  freezeout  in  the  air  flow  where  the  size-jumps  occur.  The  10 
and  100  micron  droplets  remain  completely  liquid  at  the  test  section. 
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Figure  20,  which  shows  the  predicted  water  particle  or  droplet  temperatures,  confirms 
that  the  15-80  micron  particles  have  begun  to  freeze  in  the  flow.  The  (initially)  10  micron 
particle  does  not  experience  freezing,  nor  does  the  (initially)  100  micron  particle,  up  to  their 
arrival  at  the  test  section  and,  they  have  very  different  temperatures.  The  other  particles  have 
either  become  ice  particles  (15, 20,  30  micron  particles)  or  are  slushball  or  mixture  particles  (40, 
50,  60,  80  micron  particles)  at  the  test  section.  The  results  of  the  calculation  made  with  the 
MHNT  freezing  model  indicate  the  potentially  complex  results  that  could  be  obtained  in  actual 
icing  or  adverse  weather  test  simulations  when  either  inlet  flow,  water  spray  nozzle,  or  duct 
geometry  result  in  thermal  and  kinetic  non-equilibrium  flow  conditions  at  the  test  station  or  test 

article. 

Concluding  the  results  of  test  case  3,  Figure  21  shows  the  predicted  water  particle 
velocities.  As  Figure  21  shows,  for  this  case,  all  of  the  particles,  even  the  100  micron  droplet, 
are  in,  or  very  near  to,  kinetic  or  velocity  equilibrium  with  the  air  flow  at  the  test  section  (x  =  46 

feet). 

Briefly  summarizing  this  section,  the  AEDC1DMP  code  is  capable  of  predicting  results 
of  two-phase,  dilute,  air  and  entrained  water  droplet  flows  in  research  scale  and  full-scale  icing 
and  adverse  weather  test  facilities.  The  incorporation  of  a  second  generation  or  second  order  (in 
the  hierarchy  of  physics  of  freezing  modeling)  modified  homogeneous  nucleation  theory 
(MHNT)  for  particle  freezing  in  AEDC1DMP  gives  the  code  the  capability  of  accounting  for 
particulate  freezing  over  a  wide  range  of  flow  conditions.  Appendices  C-E  describe  the  data 
input  audits  format  in  detail,  for  user  convenience.  Included  are  the  input  data  for  test  case  3. 
Appendix  F  lists  a  sample  output  from  AEDC 1 DMP  for  test  case  3 . 
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6.0  DISCUSSION  OF  RESULTS  OF  STUDY,  SUMMARY,  CONCLUSIONS,  AND 
RECOMMENDATIONS 

6. 1  Summary  of  the  Present  Study 

The  study  and  the  report  had  three  parts.  The  first  part  was  an  examination  of 
homogeneous  nucleation  theory  and  its  modifications  and  applications  to  the  prediction  of  the 
freezing  of  supercooled  pure  water  particles  and  heterogeneous  freezing  of  water  particles 
containing  motes  or  foreign  substances.  The  second  part  addressed  the  behavior  of  the 
homogenous  nucleation  rate  function,  J(T),  as  a  function  of  water  temperature.  This  function 
was  then  calibrated  to  a  set  of  heterogeneous  particle  freezing  data,  provided  by  Dorsch  and 
Hacker  [10],  by  assuming  that  the  surface  specific  free  energy  for  an  ice-water  interface  could  be 
represented  as  a  function  of  the  volume  or  diameter  of  the  particles  containing  the  motes.  With 
this  functional  form  for  the  surface  free  energy  in  hand,  the  function  J(T)  became,  in  essence, 
calibrated  to  the  Dorsch-Hacker  data  set. 

The  third  part  of  the  study  comprised  some  numerical  computations  of  two-phase,  dilute, 
air  and  entrained  water  particle  flows,  using  a  new  version  of  the  AEDC  one-dimensional, 
multiphase  flow  code,  AEDC1DMP.  The  results  of  these  calculations  demonstrated  that 
modified  homogeneous  nucleation  theory  [MHNT]  could  be  satisfactorily  used  in  numerical 
flow  models  to  identify  the  initiation  of  heterogeneous  freezing  events  in  submillimeter, 
supercooled  water  particles  in  the  computed  flows.  The  freezing  model  described  in  the  report 
accounts  for  the  effects  of  time,  particle  size,  and  particle  temperature  on  the  initiation  of 
freezing. 

6.2  Conclusions 

The  conclusions  of  the  research  reported  herein  are  as  follows: 
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(a)  Classical  homogeneous  nucleation  theory  can  be  easily  modified  to  predict 
heterogeneous  freezing  of  submillimeter,  supercooled  water  particles. 

(b)  Both  empirical  and  theoretical  methods  can  be  used  to  account  for  the  effect  of 
foreign  particles,  surface  wetting  and  dissolved  chemicals  in  the  initiation  and 
formation  of  freezing  events. 

(c)  A  particulate  freezing  model  for  heterogeneous  freezing  of  submillimeter, 
supercooled  water  particles  based  on  modified  homogeneous  nucleation  theory  has 
been  incorporated  in  a  one-dimensional,  multiphase  flow  code  (AEDC1DMP)  with 
excellent  results,  and  can  easily  be  incorporated  into  multidimensional,  multiphase 
flow  codes. 

6.3  Recommendations 

The  recommendations  of  the  present  study  are  as  follows: 

(d)  A  multidimensional,  multiphase  flow  analysis  of  water  spray  nozzle  discharge 
plumes  should  be  pursued  with  the  highest  priority.  The  analysis  should  include 
water  particle  freezing  in  the  spray  plumes  to  investigate  the  possibility  of  creating 
ice  crystals  in  the  spray  plumes.  This  topic  has  been  raised  in  the  first  part  of  this 
research  program  [1,2],  and  it  is  also  seen  as  an  important  continuation  of  the  present 
study. 

(e)  The  physics  of  water  particle  freezing  should  continue  to  be  studied  under 
heterogeneous  freezing  conditions.  The  further  development  of  models  for  the 
specific  surface  free  energy  of  the  molecular  clusters  that  trigger  freezing  events 
should  be  pursued.  Models  of  cSL,  such  as  that  described  by  Volmer  [50],  Fletcher 
[24],  Stoyanova  [36],  and  others,  should  be  developed  so  that  the  effects  on  particle 


freezing  of  different  kinds  of  motes  in  the  atmospheric  water  particles,  as  well  as  in 
ground  test  facility  water  supplies,  can  begin  to  be  addressed. 

(f)  In  view  of  the  previous  recommendation,  it  is  further  recommended  that  additional 
atmospheric  data  be  obtained  on  the  composition  of  water  particles  in  nominal  icing 
conditions.  The  data  should  include  the  chemical  and  physical  composition  of 
supercooled  water  particles  or  water  samples.  This  includes  dissolved  and/or 
entrained  chemicals,  as  well  as  the  kinds,  shapes,  and  numbers  of  solid  particles  in 
the  atmospheric  icing  environment.  Similar  effort  should  be  made  to  catalogue  the 
same  information  for  water  used  in  ground  test  simulations  of  icing  and  adverse 
weather  simulations  for  aircraft  engines  and/or  their  components. 

(g)  Consideration  of  the  freezing  process  for  supercooled  water  should  be  extended  to 
thin  films  of  water  on  aircraft  surfaces,  including  those  on  ice  accretions.  In  other 
words,  it  should  be  investigated  whether  there  are  water  run-back  modes  wherein 
thin  films  of  supercooled  water  can  be  formed  in  the  ice  accretion  process.  If  so, 
then  the  initiation  of  freezeout  of  run  back  water,  under  such  conditions,  might  be 
predicted  by  a  modified  homogeneous  nucleation  theory. 
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PLOT  OF  THE  NUMBER  OF  CRITICAL-SIZED  MOLECULAR  OUSTERS  PRESENT  IN  A  1  MICRCN  CROP 
FOR  t  =0.6  SECONDS,  AT  TEMPERATURE,  T  ,  UNDERGOING  HOMOGENEOUS  FREEZING 
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MASS  AND  ENERGY  BALANCES  FOR  TWO-PHASE  PARTICLE 

DURING  FREEZING 

B.l  The  particle  initially  starts  out  as  a  two-phase  mixture  of  water  and  ice  in 
equilibrium.  At  the  end  of  an  integration  step  (or  time,  St),  the  particle  is  still  two-phase  with  a 
different  ratio  of  ice  mass  to  liquid  mass  because  of  freezing  and  evaporation  that  occurred 
during  the  period,  St.  Heat  transfer  and  mass  transfer  effects  have  changed  the  particle  phase 
composition.  Let  the  subscripts  “I”  refer  to  the  solid  or  ice  phase,  and  “L”  refer  to  the  liquid 
phase. 

B.l. 2  Particle  Energy  Change 
The  initial  particle  energy,  E,  at  time  t,  is 

E  =  M  ,h,  +M  LhL  (1) 

The  final  particle  energy  at  time  t  +  S  Ms 

E'  =  M,'h'  +M'h'  (2) 

The  energy  change  is  given  by 

AE  =  E'-E  =  M,'h'  +ML'hl'  -M,h,+MLhL  (3) 

Because  the  process  is  isothermal, 

hj  =  h,  (4) 

hL  =hL  (5) 

Hence, 

A E  =  ^M,' -MLyL  (6) 

or 


AE  =  AM  ,h,  +  AM  LhL 


(7) 


B.  1 .3  Liquid  Mass  Change  by  Evaporation 
Evaporation  of  liquid  mass  occurs  during  the  period,  St,  in  the  amount 

SM^  =m"  AS  t 


(8) 


e 


where  A  is  the  surface  area  of  the  particle. 

For  a  moving  particle 

St  =  —  (9) 

V 

where  V  is  the  particle  velocity.  Hence, 

SMl  —  (10) 

L‘  V 


The  mass  flux  of  evaporation  is  given  by 


where 


M,  =  h 

€  m 


i-y. 


(id 


and 


K  = 


c  Lm 

p/  ' 


(12) 


h  —  average  convective  heat  transfer  coefficient  for  the  particle 
Lr  =  the  Lewis  number  for  the  particle  “film”  conditions 
cpi  =  the  specific  heat  of  the  “film”  fluid 

Y  =  mass  fraction  of  vapor  phase 

Y,  =  mass  fraction  at  particle-gas  “film”  interface 

Yco  =  free  stream  mass  fraction  of  vapor. 


Thus, 


SMLi  = 


h  Adx 


c  lmV 
pi  * 


fY  -Y  ^ 

to 

i-y. 


(13) 


* 


B.1.4  Total  Liquid  Mass  Change 

The  total  liquid  mass  change,  AMV  is  thus  in  two  parts:  a  part  due  to  freezing  of  the 
liquid,  with  subscript  Lr  and  a  part  due  to  evaporation,  with  subscript  Lt.  Hence, 


-AML=8MLf  +8MLt 

Clearly,  however, 

8MLf  =  AM; 

Thus, 

—  A Ml  =  AM;  +8MLi 

Solving, 

AML=-AM,-SML' 

B.1.5  Energy  Balance  Revisited 
The  energy  balance  can  now  be  written 

A E  =  AM  ,hj  +  AM  LhL 

or,  by  substitution 

AE-AMjh,  +(- AM,  -  SM^  )i L 
AE  =  AM ,  {h, -hJ-SM^ 
B.1.6  Heat  Balance  with  Freestream 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


The  energy  balance  for  the  particle  must  equal  the  energy  lost  or  gained  by  the  particle  to 
the  freestream.  Thus, 

AE  =  -qA8 t - m" A  hv  8 1  (21) 


where 


q=h(Ts-T_) 


(22) 


and 


hv  =  the  enthalpy  of  the  vapor  leaving  the  droplet. 

Using  previously  defined  parameters  in  equations  (10),  (11),  (13)  and  (21),  the  energy  flow  from 
the  particle  is  given  by 


A_  hAdx  /  \  hAdxh 

A E  = - -[TS-TJ- 


V 


c  LmV 

Pt  ‘ 


LzL 

1  -Y. 


(23) 


B .  1 .7  Energy  Equation  Equality 

The  two  energy  equations  for  the  particle  can  thus  be  combined  by  setting  equation  (20) 
equal  to  equation  (21) 

AE  =  £sM,(hI-hL)-8MLthL=-q  A8t-M'  Ahv8t  (24) 

Using  equation  (8)  in  equation  (24) 

AE  -  AM j  (h,  -hL)-8MLhL  =-q  A8t-8 M ^  (25) 

Rearranging  this  equation 

AM ,(h,  - hL)  =  -q  A8 1 - 8 M 4(fcv  - hL)  (26) 

with 

hv-hL=hfg  (27) 


where 


and 


M,'  =  M,+AM, 


thus, 

m'  =Mn -m'  =Ml  -8MLi -AM, 
The  mass  fraction  of  the  liquid  phase  is  then  given  by 


B.3  Particle  Properties 

The  particle  properties  for  the  two-phase  mixture  particle  of  ice  and  liquid 
given  in  terms  of  a  as  follows: 


average  density 

p= _ P±2± _ 

ap,+{l-a)pL 

average  specific  heat 

cp  =  cc  cPL  +  (l —  ot)cPI 

average  temperature 

T  =  Tl  -  T,  =  492°  R  (isothermal  assumption) 


(34) 

(35) 

(36) 

water  are 

(37) 

(38) 

(39) 


APPENDIX  C 


FORTRAN  DATA  INPUT  FORMATS  (“CARD  IMAGES”) 

The  following  data  input  statements  identify  the  data  parameters  that  are  input,  their  input 
order,  and  the  (“card  image”)  format  that  is  required  for  their  input.  The  term  LCR  is  an  acronym 
from  the  “old  days”  meaning  line  card  reader.  It  is  taken  care  of  in  the  code  compiling  process, 
and  is  of  no  concern  to  users  of  the  code. 

READ(LCR,7)  NXY,NUNITS,NCASES 
READ(LCR,707)NFREZ 
READ(LCR,8)  (XS(I),YS(I),I=1,NXY) 

READ(LCR,9)  NSTA,XPRINT,DXO 
RE  AD  (LCR,  10)  (ALP(I),I=1,10) 

NS1=NSTA+1 

READ(LCR,12)  (STA(I),I=1,NS1) 

READ(LCR,1 1)  CV,VG,TG,P 
READ(LCR,12)  (FLa),I=l,NSTA) 

READ(LCR,12)  (VL(I),I=1,NSTA) 

READ(LCR,12)  (TS(I),I=1,NSTA) 

READ(LCR,12)  (DMICRON(I),I=l,NSTA) 

The  actual  input  formats  for  the  data  are  given  below.  They  are  extremely  important 
because  the  code  AEDC1DMP  expects  to  read  in  the  data  in  the  defined  format.  Future  releases 
of  AEDC1DMP  will  use  user  friendly  data  input  formats,  but  at  present,  this  set  of  input  formats 
must  be  used. 


7  FORMAT(3I2) 

707  F0RMAT(1I2) 

8  FORMAT(2E12.0) 

9  FORMAT  ( 1 12,2E1 2.0) 

10  FORMAT(10A7) 

11  FORMAT  (4E  10.0) 

12  FORMAT(8El 0.0) 


C-l 


APPENDIX  D 

SAMPLE  DATA  INPUT  CASE 


The  data  below  are  for  a  single  data  input  case.  Multiple  cases  can  be  run  on 
AEDC1DMP  by  starting  each  new  data  set  with  card  NSTA,XPRINTJDXO  and  those  that  follow, 
through  (DMICRON(I),  1=1,  NSTA) ;  no  spaces  between  data  input  sets. 


800001 

01 

.0000 

1.5788 

.0250 

1.42511 

.0500 

1.29351 

.0750 

1.18135 

.1000 

1.08625 

.1250 

1.00599 

.1500 
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0.796079 

.2500 

0.763576 

.2750 

0.736492 

.3000 

0.713789 

.3250 

0.694557 

.3500 

0.678003 

.3750 

0.663443 

.4000 

0.650297 

.4250 

0.638074 

.4500 

0.62637 

.4750 

0.614858 

.5000 

0.603283 

.5250 

0.59145 

.5500 

0.579224 

.5750 

0.566517 

.6000 

0.55329 
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0.539538 

.6500 

0.525291 

.6750 

0.510605 

.7000 

0.495562 
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0.48026 

.7500 

0.46481 

.7750 

0.449333 

.8000 

0.433958 

.8250 
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0.38973 

.9000 
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.9500 
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.9750 
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1.0750 
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1.1000 
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1.1500 
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1.1750 
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2.0000 
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2.2000 
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2.4000 

0.283541 

2.6000 

0.283541 

2.8000 

0.283541 

100.025 

GENERIC  NOZZLE  GEOMETRY  FOR  PARAMETRIC  FREEZING  STUDY 

0.0  0.02  0.04  0.06  0.08  0.10  0.12  0.14 

0.16  0.18  2.70 

0.0002857  18.0  460.0  2116.0 

6.85E-05  2.06E-04  4.80E-04  3.62E-04  2.06E-04  3.43E-05  1.37E-05  1.37E-06 
1.37E-07  1.37E-08 


46. 

46. 

46. 

46. 

46. 

46. 

46.  46.0  46. 

46. 

543. 
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543.  543.  543. 

543. 
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APPENDIX  E 


DEFINITIONS  OF  THE  DATA  INPUT  PARAMETERS 

The  definitions  of  the  data  input  parameters  for  AEDC1DMP  are  the  same,  for  the  most 
part,  as  those  in  the  original  version  of  the  code  reported  in  Reference  53.  Their  specific 
definitions,  and  units,  where  necessary,  are  provided  below. 


NXY=  THE  NUMBER  OF  DUCT  WALL  COORDINATES  INPUT 

NUNITS=0  OR1 .  IF  NUNITS=0,  ALL  X  AND  Y  COORDINATES,  XPRINT,  DXO,  AND 
(STA(I),  1=1, NS1)  ARE  INPUT  IN  UNITS  OF  FEET.  IF  NUNITS  IS  INPUT  AS  1,  THEN  ALL 
SPACE  COORDINATES  ARE  INPUT  IN  UNITS  OF  INCHES. 

NCASES=  THE  NUMBER  OF  CASES  INPUT  TO  BE  COMPUTED. 

NFREZ=  0  OR  1 .  IF  NFREZ=0,  THE  WATER  PARTICLE  FREEZING  MODEL  IS  NOT 
ACTIVATED;  ALTERNATELY,  FOR  NFREZ=1,  THE  MODEL  WILL  BE  UTILIZED  IN 
THE  DUCT  FLOW  COMPUTATIONS. 

(XS(I),  I=1,NXY)=  THE  AXIAL  COORDINATES  OF  THE  DUCT  WALL  CURVE. 

(YS(I),  I=1,NXY)=THE  CORRESPONDING  RADIAL  COORDINATES  OF  THE  DUCT 
WALL  CURVE. 

NSTA=  THE  NUMBER  OF  WATER  INJECTION  STATIONS,  WITH  A  MAXIMUM  OF 
TEN(IO). 

XPRINT=  THE  DISTANCE  BETWEEN  AXIAL  LOCATIONS  WHERE  THE  FLOW 
SOLUTION  IS  PRINTED  OUT,  AND,  CONSEQUENTLY,  STORED  FOR  INPUT  TO  DATA 
PLOTTING  PROGRAMS  SUCH  AS  WAVEMETRICS(R)'  IGOR™ 
PRO(http://www.wavemetrics.com/). 

DXO=  THE  INITIAL  NUMERICAL  INTEGRATION  STEP  SIZE  THAT  THE  USER 
DESIRES.  BY  SETTING  DX0=0  OR  BLANK,  THE  CODE  SELECTS  A  DEFAULT  STEP 
SIZE  BASED  ON  THE  DUCT  LENGTH. 

ALP=  AN  80  COLUMN  ALPANUMERIC  FILE  IDENTIFIER. 

(STA(I),  1=1,  NS1)  =  THE  AXIAL  LOCATIONS,  OR  STATIONS,  WHERE  WATER  IS 
INJECTED  INTO  THE  FLOW,  PLUS,  THE  LAST  VALUE,  STA(NS  1),  WHICH  INDICATES 
THE  TERMINAL  AXIAL  LOCATION  FOR  FLOW  SOLUTION.  STA(NS1)  USUALLY 
EQUALS  XS(NXY)  OR  SLIGHTLY  LESS. 

CV=  FLOW  INITIAL  SPECIFIC  HUMIDITY,  LBM  OF  WATER  VAP OR/LB M  OF  DRY  AIR. 
VG=  THE  FLOW  INITIAL  VELOCITY,  FT/SECOND. 


TG=  THE  FLOW  INITIAL  STATIC  TEMPERATURE,  DEGREES  RANKINE. 

P=  THE  FLOW  INITIAL  STATIC  PRESSURE,  LBFORCE/SQUARE  FOOT. 

(FL(I),  1=1,  NSTA)  =  THE  WATER  LOADING  INJECTED  AT  EACH  STATION,  LBM  OF 
WATER/  LBM  OF  DRY  AIR. 

(VL(I),  1=1, NSTA)  =  THE  VELOCITIES  OF  THE  WATER  PARTICLES  INJECTED  AT  THE 
WATER  INJECTION  STATIONS,  FEET/SECOND. 

(TS(I),  1=1, NSTA)  =  THE  TEMPERATURES  OF  THE  WATER  PARTICLES  INJECTED  AT 
EACH  STATION,  DEGREES  RANKINE. 

(DMICRONS(I),  1=1, NSTA)  =  THE  DIAMETERS  OF  THE  WATER  PARTICLES  INJECTED 
AT  EACH  STATION,  IN  MICRONS  (1  METER= 1,000, 000  MICRONS) 

AS  A  POINT  OF  CLARIFICATION,  FL(I),VL(I),TS(I),  AND  DMICRON(I)  ALL 
CORRESPOND  TO  THE  WATER  PARTICLE  CONDITIONS  AT  THE  AXIAL  INJECTION 
STATION  LOCATED  AT  X=STA(I). 
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SAMPLE  OUTPUT  FROM  CASE  3 
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